©  Copyright  by  Jennifer  Ann  Hommema,  1999 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

j^ggg  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1999  to  00-00-1999 

4.  TITLE  AND  SUBTITLE 

Magnetomechanical  Behavior  of  Terfenol-D  Particulate  Composites 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Illinois  at  Urbana-Champaign, Theoretical  and  Applied 
Mechanics, 405  N.  Mathews  M/C  251  ,Urbana,IL, 61801 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

156 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


MAGNETOMECHANICAL  BEHAVIOR  OF  TERFENOL-D 
PARTICULATE  COMPOSITES 


BY 

JENNIFER  ANN  HOMMEMA 

B.S.,  University  of  Illinois,  1997 


THESIS 

Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Theoretical  and  Applied  Mechanics 

in  the  Graduate  College  of  the 
University  of  Illinois  at  Urbana-Champaign,  1999 


Urbana,  Illinois 


Acknowledgments 


The  author  would  like  to  thank  USACERL  for  the  funding  to  support  this  research.  Thanks 
also  goes  to  the  National  Science  Foundation  for  a  Graduate  Research  Fellowship  and  the 
Department  of  Theoretical  and  Applied  Mechanics  for  financial  support  of  the  author.  Spe¬ 
cial  thanks  goes  to  Dr.  Peter  Kurath  and  Rick  Rottet  of  AMTEL  for  aiding  the  author  in 
the  mechanical  testing  experiments.  Thanks  is  also  extended  to  Kent  Elam,  Dave  Foley,  and 
Dan  Mullis  of  the  TAM  machine  shop. 

Many  thanks  go  to  Prof.  Scott  White  of  the  Department  of  Aeronautical  and  Astro- 
nautical  Engineering  for  his  patience,  guidance,  and  encouragement  through  our  two-year 
relationship.  The  creative  and  insightful  input  of  Prof.  Nancy  Sottos  of  the  Department  of 
Theoretical  and  Applied  Mechanics  was  also  appreciated.  Prof.  James  Phillips  is  thanked 
for  the  important  role  he  has  played  in  the  author's  six-year  stay  in  the  Department  of 
Theoretical  and  Applied  Mechanics  at  the  University  of  Illinois.  He  has  greatly  aided  in  the 
development  of  the  author’s  writing  style  and  engineering  ethics,  besides  being  a  wonderfully 
caring  and  committed  mentor. 

Lastly,  the  author  wishes  to  thank  all  the  family  members  and  friends  who  listened, 
encouraged,  and  loved.  This  thesis  would  not  have  happened  without  the  support  of  these 
incredible  people. 


Table  of  Contents 


Chapter  1  Introduction  .  1 

1.1  Structural  health  monitoring  of  composite  materials .  1 

1.2  Magnetostriction .  3 

1.2.1  Ferromagnetic  “Weiss”  domain  theory  .  3 

1.2.2  Piezomagnetism .  4 

1.2.3  Direct  effect .  5 

1.2.4  Inverse  effect  .  5 

1.3  Terfenol-D .  6 

1.4  Magnetostrictive  composites .  7 

1.5  Terfenol-D  particulate  composites .  8 

1.6  Project  overview  .  9 

Chapter  2  Compression  testing .  13 

2.1  Specimen  manufacture .  13 

2.1.1  Preparation  of  Terfenol-D  .  13 

2.1.2  Terfenol-D  composite  manufacture  .  14 

2.1.3  Magnetic  flux  during  cure  and  magnetic  anneal  after  cure .  14 

2.2  Experimental  procedure  .  15 

2.3  Typical  results  and  data  reduction  .  16 

2.4  Gage  length  characterization .  17 

2.4.1  Pure  Terfenol-D .  18 

2.4.2  Terfenol-D  composites .  19 

2.4.2. 1  Terfenol-D  composites  with  chaining  in  1-direction .  20 

2. 4. 2. 2  Terfenol-D  composites  with  chaining  in  3-direction .  21 

2.4.3  Determination  of  piezomagnetic  coefficients  .  22 

2.4.3. 1  Terfenol-D .  24 

2. 4. 3. 2  Terfenol-D  composites .  24 

2.5  Effect  of  magnetic  flux  on  cyclic  response .  25 

Chapter  3  Residual  stress  analysis . 56 

3.1  Introduction .  56 

3.2  Analysis .  57 

3.3  Results .  62 

3.3.1  Material  properties .  62 


IV 


3.3.2  Thermal  mismatch .  62 

3.3.3  Chemical  shrinkage .  64 

3.3.4  Magnetostrictive  effects .  66 

3.3.5  Combination  of  all  effects  .  67 

3.4  Summary  .  68 

Chapter  4  Finite  element  analysis . 78 

4.1  Introduction .  78 

4.1.1  Analogy  between  piezomagnetism  and  piezoelectricity .  78 

4.1.2  Boundary  conditions .  80 

4.1.3  Symmetry  from  polarization .  81 

4.2  Finite  element  formulation .  83 

4.2.1  Preprocessing  with  PATRAN  .  83 

4.2.2  Analysis  with  ABAQUS .  84 

4.2.3  Modification  of  the  input  hie .  85 

4.3  Composite  cylinder  model .  87 

4.3.1  Geometry,  loading,  and  boundary  conditions .  87 

4.3.2  Material  properties .  87 

4.4  Convergence  study .  88 

4.5  Results .  89 

Chapter  5  Summary  and  conclusions . 104 

Appendix  A  Gage  length  characterization  results . 107 

Appendix  B  Cyclic  testing  results . 129 

Appendix  C  Sample  input  file  for  FEA . 137 

References . 140 


v 


List  of  Tables 


1.1  Simplification  of  tensor  notation  to  matrix  form  of  piezomagnetic  equations.  5 

2.1  Comparison  of  piezomagnetic  coefficients  for  the  different  Terfenol-D  compos¬ 
ites  discussed  in  Sections  2.4.2. 1  and  2. 4. 2. 2 .  25 

2.2  Comparison  of  Af?“ax  for  different  combinations  of  curing  and  annealing  mag¬ 
netic  fluxes  for  (1,1,1)  Terfenol-D  composites .  27 

2.3  Comparison  of  Af?“ax  for  different  combinations  of  curing  and  annealing  mag¬ 
netic  fluxes  for  (3,3,1)  Terfenol-D  composites .  27 

3.1  Summary  of  material  properties  used  in  modified  Eshelby  model  for  residual 

stresses .  62 

3.2  Summary  of  chemical  shrinkage  experiment  results .  66 

3.3  Summary  of  residual  stresses  for  spherical  particle  model  and  ellipsoidal  par¬ 
ticle  chain  model .  68 

4.1  Comparison  of  piezomagnetic  and  piezoelectric  governing  equations .  79 

4.2  Comparison  of  piezomagnetic  and  piezoelectric  quantities  with  units .  79 

4.3  Comparison  of  interface  conditions  for  piezomagnetism  and  piezoelectricity.  .  80 

4.4  Piezoelectric  elements  available  in  ABAQUS .  84 

4.5  Summary  of  material  properties  used  in  finite  element  analysis  of  composite 

cylinder  model .  88 


vi 


List  of  Figures 


1.1  Schematic  of  tagging  concept  with  magnetostrictive  particulates .  11 

1.2  Direct  and  inverse  effects  of  magnetostriction .  11 

1.3  Domain  evolution  under  applied  magnetic  held  [32] .  12 

2.1  Geometry  of  uniaxial  compression  specimen .  28 

2.2  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  with 

a  magnetic  flux  through  the  length .  29 

2.3  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  with 

a  magnetic  flux  through  the  width .  29 

2.4  Experimental  setup  for  uniaxial  compression  tests .  30 

2.5  Schematic  of  Hall  probe  used  to  measure  magnetic  flux .  30 

2.6  Typical  response  of  (1,1,1)  composite  specimen  under  uniaxial  compressive 

stress .  31 

2.7  Typical  response  of  (1,1,1)  composite  specimen  under  uniaxial  compressive 

stress  with  magnetic  flux  zeroed  by  subtracting  magnetic  flux  value  at  zero 
load,  Djnit .  32 

2.8  Distance  dependence  of  A B™*  measurement  on  a  (3,3,1)  Terfenol-D  com¬ 
posite  specimen  loaded  to  —10  MPa  and  annealed  at  —800  G  before  each 

test .  33 

2.9  Compression  specimen  marked  with  grid  for  gage  length  characterization.  .  .  34 

2.10  Typical  gage  length  response  of  a  (1,1,1)  Terfenol-D  composite .  35 

2.11  Typical  gage  length  response  of  a  (1,1,1)  Terfenol-D  composite  with  magnetic 

flux  zeroed  by  subtracting  .  36 

2.12  Typical  ambient  magnetic  flux  measured  on  a  polyester  specimen .  37 

2.13  Comparison  of  AHmax,  the  change  in  magnetic  flux  over  —40  MPa,  for  a  pure 
Terfenol-D  specimen  annealed  at  a  magnetic  flux  of  —800  G  in  the  1-clirection.  38 

2.14  Comparison  of  the  change  in  magnetic  flux  over  —40  MPa,  AHmax,  for  a  pure 

Terfenol-D  specimen  annealed  at  a  magnetic  flux  of  —800  G  in  the  ^-direction, 
where  z  is  the  axis  along  the  length  of  the  specimen  and  r  points  outward.  .  39 

2.15  Magnetic  flux  lines  of  a  uniformly  magnetized  cylinder  [26] .  40 

2.16  Comparison  of  AHmax,  the  change  in  magnetic  flux  over  —10  MPa,  for  a 

(1,1,1)  Terfenol-D  composite  specimen  cured  and  annealed  at  —800  G.  .  .  .  41 

2.17  Comparison  of  AHmax,  the  change  in  magnetic  flux  over  —10  MPa,  for  a  (1,3,1) 

Terfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G.  .  42 

vii 


2.18  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  (1,1,1) 

and  (1,3,1)  Terfenol-D  composite  specimens  cured  and  annealed  at  a  magnetic 
flux  of  —800  G,  where  z  is  the  axis  along  the  length  of  the  specimen  and  r 
points  outward .  43 

2.19  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  a  (3,1,1) 

Terfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G.  .  44 

2.20  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  a  (3,2,1) 

Terfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G.  .  45 

2.21  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  a  (3,3,1) 

Terfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G.  .  46 

2.22  Geometries  of  compression  specimens  converted  to  coordinate  system  used  in 

piezomagnetic  literature .  47 

2.23  Piezomagnetic  model  prediction  and  experimental  data  for  Terfenol-D  speci¬ 
men  poled  in  the  same  direction  as  the  compressive  stress .  48 

2.24  Piezomagnetic  model  prediction  and  experimental  data  for  (1,1,1)  Terfenol-D 

composite  specimen .  49 

2.25  Piezomagnetic  model  prediction  and  experimental  data  for  (1,3,1)  Terfenol-D 

composite  specimen .  50 

2.26  Piezomagnetic  model  prediction  and  experimental  data  for  (3,1,1)  Terfenol-D 

composite  specimen .  51 

2.27  Piezomagnetic  model  prediction  and  experimental  data  for  (3,2,1)  Terfenol-D 

composite  specimen .  52 

2.28  Piezomagnetic  model  prediction  and  experimental  data  for  (3,3,1)  Terfenol-D 

composite  specimen .  53 

2.29  Typical  cyclic  response  for  a  Terfenol-D  composite  loaded  in  uniaxial  com¬ 
pression  to  —10  MPa  with  annealing  only  once  before  the  first  cycle .  54 

2.30  Typical  cyclic  response  for  a  Terfenol-D  composite  loaded  in  uniaxial  com¬ 
pression  to  —10  MPa  with  annealing  only  once  before  the  first  cycle .  55 

3.1  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  un¬ 
der  a  magnetic  field  in  X%  direction  [2] .  70 

3.2  Geometry  of  the  ellipsoidal  inclusion  in  the  modified  Eshelby  model .  71 

3.3  Model  geometries  for  the  modified  Eshelby  model .  72 

3.4  Residual  stresses  from  thermal  mismatch .  73 

3.5  Temperature  recorded  by  thermocouple  at  center  of  polyester  specimen  cured 

at  room  temperature  between  permanent  magnets  at  800  G .  74 

3.6  Residual  stresses  from  chemical  shrinkage .  75 

3.7  Residual  stresses  from  magnetostrictive  strains .  76 

3.8  Direct  effect  of  magnetostriction  for  a  Terfenol-D  rod  with  no  prestress  [21].  77 

4.1  Geometry  of  the  composite  cylinder  model .  91 

4.2  Mesh  1  with  14,  869  nodes  and  4,  860  elements .  92 

4.3  Mesh  2  with  17,  750  nodes  and  5,  813  elements .  93 

4.4  Difference  in  magnetic  flux,  Bi,  across  interface  at  lower  10%  of  the  composite 

cylinder  model .  94 

viii 


4.5  Difference  in  magnetic  flux,  B i,  across  interface  at  upper  10%  of  the  composite 

cylinder  model .  94 

4.6  Difference  in  shear  stress,  a\2)  across  interface  at  lower  10%  of  the  composite 

cylinder  model .  95 

4.7  Difference  in  shear  stress,  <j12,  across  interface  at  upper  10%  of  the  composite 

cylinder  model .  95 

4.8  Contour  plot  of  the  magnetic  flux,  B i,  in  the  Terfenol-D  rod  of  the  composite 

cylinder  model .  96 

4.9  Contour  plot  of  the  magnetic  flux,  B2,  in  the  Terfenol-D  rod  of  the  composite 

cylinder  model .  97 

4.10  Contour  plot  of  the  magnetic  flux,  Bi,  in  the  polyester  matrix  of  the  composite 

cylinder  model .  98 

4.11  Contour  plot  of  the  magnetic  flux,  B2,  in  the  polyester  matrix  of  the  composite 

cylinder  model .  99 

4.12  The  magnetic  flux,  B\ ,  along  the  line  y  =  5  mm  of  the  composite  cylinder 

model .  100 

4.13  The  magnetic  flux,  B2:  along  the  line  y  =  5  mm  of  the  composite  cylinder 

model .  101 

4.14  Magnetic  flux,  Bf ,  along  the  interface  in  the  Terfenol-D  rod .  102 

4.15  Magnetic  flux,  B J,  along  the  interface  in  the  Terfenol-D  rod .  103 

A.l  Magnetic  flux  ,B i,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 

flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position . 108 

A. 2  Zeroed  magnetic  flux,  Bl}  over  the  gage  length  of  pure  Terfenol-D  specimen, 

annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position.  108 
A. 3  Magnetic  flux,  B2 ,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 

flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position . 109 

A. 4  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  pure  Terfenol-D  specimen, 

annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position.  109 
A. 5  Magnetic  flux,  B2,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 

flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position . 110 

A. 6  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  pure  Terfenol-D  specimen, 

annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position.  110 
A. 7  Magnetic  flux,  B%,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 

flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position . Ill 

A. 8  Zeroed  magnetic  flux,  B%,  over  the  gage  length  of  pure  Terfenol-D  specimen, 

annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position.  Ill 
A. 9  Magnetic  flux,  B\ ,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  112 

A. 10  Zeroed  magnetic  flux,  B i,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position .  112 


IX 


A.  11  Magnetic  flux,  B2,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  113 

A. 12  Zeroed  magnetic  flux,  B2)  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  113 

A. 13  Magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,1,1)  Terfenol-D 
composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G 

applied  before  testing  at  each  position .  114 

A.  14  Zeroed  magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,1,1) 
Terfenol-D  composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of 

—800  G  applied  before  testing  at  each  position . 114 

A.  15  Magnetic  flux,  H3,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  115 

A. 16  Zeroed  magnetic  flux,  B$,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  115 

A. 17  Magnetic  flux,  Bi,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 

imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  116 

A. 18  Zeroed  magnetic  flux,  B 1,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  116 

A.  19  Magnetic  flux,  B2,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 

imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  117 

A. 20  Zeroed  magnetic  flux,  B2 ,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  117 

A. 21  Magnetic  flux,  B2 ,  over  the  gage  length  on  the  side  face  of  (1,3,1)  Terfenol-D 

composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G 

applied  before  testing  at  each  position .  118 

A. 22  Zeroed  magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,3,1) 
Terfenol-D  composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of 

—800  G  applied  before  testing  at  each  position . 118 

A. 23  Magnetic  flux,  B$,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  119 

A. 24  Zeroed  magnetic  flux,  H3,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position .  119 


x 


A. 25  Magnetic  flux,  B i,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  120 

A. 26  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  120 

A. 27  Magnetic  flux,  B2,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  121 

A. 28  Zeroed  magnetic  flux,  B2)  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  121 

A. 29  Magnetic  flux,  B3,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  122 

A. 30  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  122 

A. 31  Magnetic  flux,  B{l  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 

imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  123 

A. 32  Zeroed  magnetic  flux,  B i,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  123 

A. 33  Magnetic  flux,  B2,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  124 

A. 34  Zeroed  magnetic  flux,  B2)  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  124 

A. 35  Magnetic  flux,  B3,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  125 

A. 36  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  125 

A. 37  Magnetic  flux,  B{,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 

imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  126 

A. 38  Zeroed  magnetic  flux,  B\,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position .  126 


xi 


A. 39  Magnetic  flux,  B2,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  127 

A. 40  Zeroed  magnetic  flux,  B2)  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 

testing  at  each  position .  127 

A. 41  Magnetic  flux,  B5,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 

imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position .  128 

A.  42  Zeroed  magnetic  flux,  B%,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 

specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position .  128 

B. l  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 

ing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  130 

B.2  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  130 

B.3  Cyclic  magnetic  response,  B\,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —800  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 

only .  131 

B.4  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —800  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 

only .  131 

B.5  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  132 

B.6  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  132 

B.7  Cyclic  magnetic  response,  B±,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 

only .  133 

B.8  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 

only .  133 

B.9  Cyclic  magnetic  response,  B:> ,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  134 

B.10  Cyclic  magnetic  response,  B:> ,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 
only .  134 


xii 


B.ll  Cyclic  magnetic  response,  B3,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  135 

B.12  Cyclic  magnetic  response,  B3,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle 

only .  135 

B.13  Cyclic  magnetic  response,  /i 3 ,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 

only .  136 

B.14  Cyclic  magnetic  response,  B3,  of  (3,3,1)  Terfenol-D  composite  specimen,  cur¬ 
ing  flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle 
only .  136 


xiii 


Chapter  1 
Introduction 


1.1  Structural  health  monitoring  of  composite 
materials 

A  structural  composite  is  defined  as  a  material  consisting  of  two  or  more  phases  on  a  macro¬ 
scopic  scale.  The  properties  of  the  composite  material  are  improved  from  the  properties  of 
the  constituent  phases  acting  alone.  Typically,  a  structural  composite  consists  of  a  rein¬ 
forcement  phase,  a  matrix  phase,  and  an  interphase  or  interface.  The  reinforcement  phase  is 
usually  stiffer  and  stronger  and  takes  the  form  of  particulates  or  fibers.  The  matrix  could  be 
a  ceramic,  a  metal,  or  a  polymer.  One  of  the  first  documented  uses  of  structural  composites 
is  straw-reinforced  clay  bricks  in  ancient  Egypt  [6]*. 

More  recently,  in  the  twentieth  century,  structural  polymer  composites  have  found  many 
applications  in  the  aerospace,  automotive,  sporting  goods,  and  medical  industries.  Struc¬ 
tural  polymer  composites  have  high  specific  stiffness  (modulus  to  density  ratio)  and  high 
specific  strength  (strength  to  density  ratio).  Some  polymer  composites  also  have  increased 
corrosion  resistance,  wear  resistance,  thermal  stability,  and/or  fatigue  life.  Advanced  manu¬ 
facturing  methods  such  as  autoclave  molding,  filament  winding,  pultrusion,  and  resin  transfer 
molding  allow  structural  composites  to  be  designed  with  complex  geometries.  One  of  the 
*References  are  listed  alphabetically  by  author  at  the  end  of  the  thesis,  beginning  on  page  140. 
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disadvantages  of  polymer  composites  is  their  failure  process.  The  most  common  matrices 
are  thermosetting  polymers  such  as  epoxies  and  polyesters,  which  fail  in  a  brittle  manner. 
Brittle  failure  of  structural  polymer  composites  is  often  catastrophic  and  occurs  with  little 
warning.  For  example,  when  a  multidirectional  laminate  experiences  impact  loading,  delam¬ 
inations  and  other  damage  occur  beneath  the  surface  of  the  composite.  This  damage  can 
not  be  detected  by  visual  inspection  alone  and  could  lead  to  unexpected  failure  [6]. 

New  techniques  are  necessary  to  monitor  the  health  of  a  structure.  These  techniques 
should  detect  damage  and  prevent  failure,  or  at  least  anticipate  it.  There  are  many  different 
ways  to  achieve  structural  health  monitoring  through  nondestructive  evaluation  including 
visual  inspection,  ultrasonic  C-scan,  and  x-ray  radiography.  A  recently  developed  technique 
of  interest  for  polymer  composites  is  tagging,  in  which  a  smart  material  phase  is  added  to 
the  matrix  of  the  composite  as  shown  in  Figure  1.1.  The  presence  of  the  smart  material 
allows  the  composite  to  be  self-assessing.  A  self-assessing  material  is  able  to  interrogate  and 
assess  its  own  state  of  health.  The  tagging  phase  could  take  the  form  of  thin  films,  wires, 
or  particles.  Common  classes  of  smart  materials  which  could  be  used  include  piezoelectric, 
magnetostrictive,  and  shape-memory  alloy. 

Magnetostrictive  tagging  offers  unique  capabilities  for  self-assessment.  Utilizing  the 
so-called  inverse  effect,  a  magnetostrictive  material  develops  or  changes  its  magnetic  held 
under  the  application  of  stresses.  Thus,  by  measuring  the  magnetic  fields  surrounding  a 
magnetostrictive  material,  a  measure  of  the  stresses  on  that  material  is  obtained.  If  the 
stresses  are  elevated  due  to  the  presence  of  damage,  the  magnetic  held  will  also  be  elevated. 
Self-assessment  is  achieved  by  incorporating  magnetostrictive  tags  into  a  structural  material 
and  monitoring  the  magnetic  held  emanating  from  the  material. 

This  thesis  investigates  the  magnetomechanical  behavior  of  Terfenol-D  particulate  com¬ 
posites.  Terfenol-D  is  a  magnetostrictive  material.  By  grinding  the  material  into  small 
particles,  the  powder  is  easily  mixed  into  a  polymer  resin  during  manufacture  of  the  com¬ 
posite.  Once  the  resin  is  cured,  the  composite  material  has  self-assessing  capabilities  since 
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the  magnetostrictive  particles  respond  to  changes  in  the  stress  state. 


1.2  Magnetostriction 

Magnetostriction  is  defined  as  the  deformation  of  any  substance  clue  to  the  presence  of  mag¬ 
netic  fields.  It  occurs  in  ferromagnetic  materials  including  iron,  cobalt,  and  nickel  and  was 
first  documented  by  Joule  in  1842  in  iron  bars.  The  two  types  of  magnetostriction  are  vol¬ 
ume  magnetostriction  and  Joule  magnetostriction.  Volume  magnetostriction  is  an  isotropic 
volume  effect,  while  Joule  magnetostriction  is  anisotropic.  Referring  to  Figure  1.2,  consider 
a  magnetostrictive  sphere  in  the  absence  of  any  magnetic  field  (H).  When  a  magnetic  field 
is  applied,  the  sphere  elongates  into  an  ellipsoid  with  the  symmetry  axis  along  the  direc¬ 
tion  of  the  applied  magnetic  field.  This  process  is  called  Joule  magnetostriction.  There 
is  a  limit  to  this  induced  strain,  which  is  known  as  the  saturation  magnetostriction.  The 
thought  experiment  above  demonstrates  the  direct  effect  of  magnetostriction,  but  an  inverse 
effect  also  exists.  In  the  inverse  effect,  mechanical  stresses  (T)  on  a  body  affect  its  mag¬ 
netic  state  [32],  These  two  effects  are  shown  in  Figure  1.2,  where  the  top  spheres  show  the 
mechanical  deformation  and  the  lower  spheres  show  the  magnetic  orientation,  as  represented 
by  arrows. 

1.2.1  Ferromagnetic  “Weiss”  domain  theory 

Pierre  Weiss  defined  a  domain  in  a  ferromagnetic  material  as  a  volume  in  which  all  the 
elementary  magnetic  moments  are  aligned  in  the  same  direction  [32],  The  moments  are 
aligned  with  one  of  the  “easy  axes” ,  which  are  defined  by  the  crystal  structure.  Domains  are 
separated  by  Bloch  walls  in  which  the  magnetization  rotates  from  one  easy  axis  to  another 
easy  axis.  In  a  demagnetized  body,  the  net  magnetization  is  zero  as  shown  in  (a)  of  Figure  1.3. 
The  domains  are  labeled  Dj  and  each  has  an  arrow  representing  the  direction  of  the  magnetic 
moments  or  the  magnetization.  When  a  magnetic  field  is  applied  to  the  domains,  several 
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different  domain  motions  can  occur.  First,  the  domains  rotate  180°  to  align  with  the  applied 
held  as  shown  in  (b).  Next,  90°  domain  wall  displacement  occurs.  In  each  of  these  processes, 
Bloch  walls  disappear  when  the  domains  merge.  Finally,  the  magnetization  will  rotate  to 
the  same  direction  as  the  applied  held  [32], 

1.2.2  Piezomagnetism 

In  general,  magnetostriction  is  a  second-order,  nonlinear  effect,  but  it  can  be  treated  as  a 
first-order,  linear  effect  over  a  certain  range  of  operation  [22],  The  linear  magnetostrictive 
effect  is  called  piezomagnetism  and  is  described  by  the  equations, 


Sij  T  dkijHfc  (1.1a) 

Bi  =  dmTki  +  tfkHk  (1-lb) 


where  TVJ  is  the  stress  tensor,  S%3  is  the  strain  tensor,  Bi  is  the  magnetic  flux  vector,  and  Hi 
is  the  magnetic  held  vector.  The  tensor,  dijk,  is  composed  of  material  properties  known  as 
piezomagnetic  coefficients.  The  tensor  equations  can  be  simplified  to  the  matrix  equation 
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Table  1.1:  Simplification  of  tensor  notation  to  matrix  form  of  piezomagnetic  equations. 


ij  and  kl  (tensor  indices) 

11 

22 

33 

23  and  32 

31  and  13 

12  and  21 

p  and  q  (matrix  form) 

1 

2 

3 

4 

5 

6 

1.2.3  Direct  effect 

The  direct  effect  is  very  well  documented  in  the  literature.  Much  research  has  been  conducted 
to  optimize  the  properties  of  Terfenol-D  for  magnetostrictive  transducers  and  actuators. 
Most  of  this  transducer  research  uses  one-dimensional  piezomagnetic  theory.  The  piezomag¬ 
netic  coefficient,  d 33,  is  reported  to  have  a  value  between  1.0  x  10-8  and  2.5  x  10-8  m/A 
by  Moffett  et  al.  [20].  A  similar  transducer  experiment  yielded  lower  d33  values  between 
2.5  x  1(T9  and  5.5  x  1(T9  m/A  [7]. 

1.2.4  Inverse  effect 

The  inverse  effect  of  magnetostriction  is  sometimes  referred  to  as  the  Villari  effect  [32],  There 
is  significantly  less  published  work  on  the  inverse  effect  than  the  direct  effect.  The  Villari 
effect  in  ferrites  has  been  proposed  as  a  possible  means  of  measuring  stresses  [1].  Steel  and 
nickel  samples  were  externally  stressed,  yielding  changes  in  the  magnetic  field  measured  with 
a  Permalloy  magnetoresistor  [17,  18]. 

Within  the  large  published  accounts  of  the  direct  effect,  there  is  some  discussion  of 
the  inverse  effect  in  a  very  specialized  application.  Terfenol-D  actuators  are  compressively 
prestressed  before  application  of  the  driving  magnetic  field.  This  prestress  increases  the 
saturation  magnetostriction  and  the  piezomagnetic  coefficient  in  the  direction  of  the  applied 
magnetic  field  [5].  It  is  hypothesized  that  the  compressive  prestress  rotates  the  domains  away 
from  their  initial  poling  axis  along  the  length  axis  of  the  actuator  rod.  Then,  the  application 
of  the  magnetic  field  can  rotate  more  domains  into  the  direction  of  the  field  along  the  length 
axis,  resulting  in  a  higher  magnetostrictive  strain  [32], 
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1.3  Terfenol-D 


Terfenol-D  is  a  ternary  alloy  of  terbium,  dysprosium,  and  iron  with  the  chemical  formula 
Tbo  3Dyo.7Fe2.o-  It  is  called  a  “giant  magnetostrictive” ,  since  it  is  capable  of  a  large  satura¬ 
tion  magnetostrictive  strain.  The  saturation  magnetostriction  of  Terfenol-D,  as  defined  in 
Section  1.2  and  denoted  by  As,  has  been  measured  to  have  a  value  as  high  as  2.4  x  10“3. 
This  magnetostrictive  strain  is  the  largest  known  room  temperature  saturation  magnetostric¬ 
tion  [32],  Values  for  the  saturation  magnetostriction  reported  by  the  material  manufacturer, 
Etrerna  Products,  are  between  1.5  x  10~3  and  2.0  x  10-3  [8].  For  rods  grown  by  the  Bridg¬ 
man  technique  and  the  free-standing-zone  melt  technique,  As  values  were  1.6  x  10~3  and 
2.0  x  10~3,  respectively.  For  a  single-crystal  rod,  As  was  2.4  x  10-3  [11]. 

The  reported  material  properties  of  Terfenol-D  vary  widely.  This  variation  is  in  part  due 
to  the  dependence  of  the  properties  on  the  crystal  structure  and  the  magnetic  and  stress 
conditions  under  which  the  properties  are  measured  [32],  For  three  Terfenol-D  rods  grown 
by  the  Bridgman,  the  free-standing-zone  melt,  and  the  single-crystal  techniques,  the  Young’s 
modulus  was  consistently  100  GPa  [11].  Varying  the  magnitude  of  the  magnetic  flux  on  a 
rod  from  0  to  4.5  kG  caused  the  Young’s  modulus  to  increase  from  44  GPa  to  110  GPa  [4], 
In  an  article  on  a  magnetostrictive  finite  element  analysis,  Kannan  and  Dasgupta  selected 
30  GPa  and  0.444  for  the  Young’s  modulus  and  Poisson  ratio,  respectively  [19].  Values  for 
tensile  and  compressive  strength  are  28  and  —700  MPa,  respectively  [32], 

The  piezoelectric  properties  of  Terfenol-D  have  been  measured  for  the  direct  effect  with 
a  one- dimensional  model;  therefore,  only  the  d33  value  is  found  in  the  literature.  In  a 
review  article,  reported  ^33  values  were  as  small  as  10  x  10~9  m/A  or  as  large  as  90  x 
10-9  m/A  depending  on  the  magnitude  of  the  magnetic  field  and  prestress  used  in  the 
actuator  experiment  [14],  For  a  Bridgman  technique  rod,  d33  was  24  x  10~9  m/A  and  for  a 
free-standing-zone  melt,  ^33  was  57  x  10~9  m/A  [32],  Variation  of  compressive  prestress  and 
magnetic  bias  field  showed  d33  decreasing  from  25  x  10~9  to  4  x  10~9  m/A  as  the  bias  field 
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and  prestress  were  increased  [20]. 


The  values  for  permeability  of  Terfenol-D  are  reported  as  relative  values,  where 


rel  _  /% 

Vij  ~  .. 
ho 


(1.3) 


where  /i0  is  the  permeability  of  free  space.  A  range  of  4  to  17  was  measured  for  the  relative 
permeability  as  a  function  of  frequency  [9].  Etrerna  Products  reported  relative  permeability 
values  between  5  and  10  [32],  By  variation  of  compressive  prestress  and  magnetic  bias  field, 
relative  permeability  measured  for  a  Terfenol-D  rod  was  between  4  and  15,  with  the  highest 
value  corresponding  to  the  lowest  prestress  and  lowest  magnetic  bias  field  [20]. 


1.4  Magnetostrictive  composites 

Magnetostrictive  composites  have  been  studied  and  used  in  direct  effect  applications  [30,  13]. 
Rod-shaped  particulate  composites  with  a  nonmetallic  insulating  matrix  were  manufactured 
with  70-80%  volume  fraction  of  Terfenol-D  [30].  In  a  transducer  experiment,  the  rod  was 
placed  in  a  solenoid  that  provided  a  driving  magnetic  field.  Strain  gages  were  used  to  measure 
longitudinal  strain  and  an  induction  pick-up  coil  measured  the  magnetization  changes.  One 
type  of  composite  was  isotropic,  while  the  second  type  was  anisotropic  because  the  particles 
had  been  magnetically  oriented.  The  anisotropic  rods  performed  slightly  better  than  the 
isotropic  rods.  The  saturation  magnetostrictions  measured  for  the  rods  were  approximately 
40%  of  the  saturation  magnetostriction  of  pure  Terfenol-D.  The  piezomagnetic  constant,  ^33, 
was  measured  to  be  between  3  x  10~9  and  6  x  10”9  m/A  [30]. 

Using  a  similar  experimental  procedure,  Hudson  et  al.  investigated  the  dependence  of 
magnetostrictive  composite  properties  on  particle  size  and  volume  fraction.  Composites 
consisted  of  an  epoxy  matrix  and  particles  which  ranged  in  size  from  106-710  /im.  Only 
three  volume  fractions  were  tested:  55%,  62%  and  70%.  The  first  conclusion  was  that  the 
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magnetomechanical  coupling  coefficient  was  independent  of  particle  size  and  volume  fraction. 
Elastic  moduli  were  independent  of  particle  size  but  dependent  on  volume  fraction  [13]. 

Jolly  et  al.  studied  the  magnetostrictive  direct  effect  with  an  application  to  elastomer 
composites  [15].  Carbonyl  iron  particles  were  embedded  in  a  polymer  matrix.  The  particles 
were  exposed  to  a  strong  magnetic  held  during  cure  which  caused  the  particle  to  migrate  into 
particle  chains  parallel  with  the  applied  held.  When  the  elastomer  composite  was  cured,  it 
was  tested  in  a  double  lap  shear  to  determine  how  much  the  stiffness  changed  with  variation 
of  an  applied  magnetic  held  [15].  Metal  matrix  magnetostrictive  composites  were  successfully 
fabricated  by  Pinkerton  et  al.  [25].  The  composites  consisted  of  a  matrix  of  either  iron  or 
aluminum  and  a  50%  volume  fraction  of  SmFe2  ribbons.  These  composites  exhibited  large 
saturation  magnetostrictive  strains  and  good  mechanical  integrity. 

1.5  Terfenol-D  particulate  composites 

For  the  purpose  of  structural  health  monitoring,  Terfenol-D  particulate  composites  were  hrst 
proposed  by  Rogers  et  al.  [29].  These  researchers  were  able  to  demonstrate  elevated  magnetic 
holds  at  the  location  of  a  delamination  in  a  Terfenol-D  tagged  epoxy  composite.  They  also 
investigated  the  magnetic  helds  measured  in  a  composite  consisting  only  of  Terfenol-D  and 
epoxy.  A  similar  report  was  made  by  Jones  and  Liang  [16].  Their  paper  also  presented 
results  showing  that  embedded  Terfenol-D  particles  could  be  used  to  locate  a  delamination 
in  a  composite.  Unfortunately,  no  information  was  included  as  to  particle  size,  distribution, 
or  manufacture  of  the  specimen. 

An  initial  study  was  conducted  by  White  and  Albers  to  determine  the  feasibility  of  using 
magnetostrictive  Terfenol-D  tags  for  structural  health  monitoring  [36].  The  chaining  phe¬ 
nomenon  documented  by  Jolly  et  al.  [15]  was  observed  in  the  polymer  composites  of  White 
and  Albers.  They  documented  a  definite,  repeatable,  measurable  magnetic  response  under 
tensile  loading  for  several  different  composites.  They  also  showed  that  elevated  magnetic 


signals  were  seen  at  the  site  of  a  delamination  [36]. 

This  initial  study  was  expanded  by  Brouwers  [2],  who  conducted  extensive  mechanical 
testing  on  Terfenol-D/polyester  composites  with  and  without  glass  fiber  reinforcement.  His 
experimental  work  used  a  Hall  probe  to  detect  changes  in  magnetic  flux  as  a  composite  was 
loaded  in  tension.  He  demonstrated  that  the  application  of  a  magnetic  field  during  cure 
caused  the  chaining  phenomenon  and  increased  the  magnetic  response  of  the  composites. 
He  also  introduced  the  concept  of  “magnetic  annealing”  to  increase  the  magnetic  response 
of  the  composites.  Damage  detection  tests  for  the  Terfenol-D/polyester  composites  were 
quite  successful.  Magnetic  flux  levels  were  elevated  at  the  edge  of  a  hole  and  also  at  the 
tip  of  a  square  notch.  Similar  tests  on  Terfenol-D/glass/polyester  composites  were  inconclu¬ 
sive.  Unfortunately,  only  a  single  component  of  the  magnetic  flux  was  reported  in  this  work. 
Three-dimensional  measurements  of  the  magnetic  flux  were  obtained  later  and  reported  [35] . 
The  most  significant  result  in  this  latter  work  is  that  the  largest  maximum  magnetic  response 
was  in  the  direction  of  the  curing  and  annealing  magnetic  fluxes.  These  responses  were  the 
largest  ever  reported  for  Terfenol-D  particulate  composites.  The  piezomagnetic  coefficients 
measured  were  <^33  equal  to  7.97  x  1CT11  m/A  and  d^i  equal  to  2.49  x  1CD11  m/A  [34], 
Overviews  of  these  magnetostrictive  tagging  results  as  they  apply  to  structural  health  mon¬ 
itoring  are  given  by  Quattrone  et  al.  and  White  [27,  33]. 

1.6  Project  overview 

The  research  in  this  thesis  is  divided  into  three  distinct  areas.  First,  the  magnetomechanical 
behavior  of  pure  Terfenol-D  and  Terfenol-D  composites  is  studied  in  uniaxial  compression  in 
Chapter  2.  The  Terfenol-D  composites  tested  consist  of  a  polyester  matrix  and  2.2%  volume 
fraction  Terfenol-D  particulates.  The  three-dimensional  magnetic  response  is  measured  over 
the  gage  length  of  the  specimens  in  uniaxial  compression.  Also,  the  piezomagnetic  coefficients 
are  predicted  for  the  pure  Terfenol-D  specimen  and  the  Terfenol-D  composites.  Cyclic  testing 
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is  presented,  showing  the  stability  of  the  response  after  multiple  loading  and  unloading  cycles. 
Different  combinations  of  curing  magnetic  flux  and  annealing  magnetic  flux  are  compared. 
In  Chapter  3,  a  modified  Eshelby  model  is  used  to  predict  the  residual  stresses  that  develop 
during  the  manufacturing  process  of  the  Terfenol-D  composites.  Finally,  a  finite  element 
analysis  of  a  Terfenol-D  rod  embedded  in  a  polyester  cylinder  is  presented  in  Chapter  4. 
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Figure  1.1:  Schematic  of  tagging  concept  with  magnetostrictive  particulates. 


Figure  1.2:  Direct  and  inverse  effects  of  magnetostriction. 
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Figure  1.3:  Domain  evolution  under  applied  magnetic  field  [32], 


12 


Chapter  2 

Compression  testing 


2.1  Specimen  manufacture 

2.1.1  Preparation  of  Terfenol-D 

One  pure  Terfenol-D  specimen  (12.5  x  12.5  x  25  mm)  was  purchased  from  Etrema  Products, 
Inc.  (Ames,  1A).  The  specimen  was  machined  from  a  near-single  crystal  rod  with  its  growth 
axis  in  the  1-direction  as  shown  in  Figure  2.1. 

To  produce  the  Terfenol-D  powder,  amorphous  Terfenol-D  fingers  were  purchased  from 
Etrema  Products,  Inc.  (Ames,  IA).  These  fingers  were  ground  in  an  argon  environment 
and  screened  to  30  mesh.  Once  the  particles  were  all  30  mesh  or  smaller,  they  were  stored 
under  argon  and  sent  to  Union  Process,  Inc.  (Akron,  OH).  At  Union  Process,  the  particles 
were  ball  milled  under  argon.  Machine  cutting  oil  was  used  to  prevent  the  particles  from 
sticking  together.  Particles  that  were  milled  without  oil  clumped  together  and  exhibited  less 
magnetostrictive  behavior  when  tested  compared  with  particles  milled  under  oil;  therefore 
only  particles  milled  in  oil  were  used  for  the  composites  in  this  study.  The  particles  in  oil 
were  shipped  back  in  an  argon  environment.  Vacuum  filtration  with  a  silica  frit  was  used  to 
separate  the  powder  from  the  oil.  The  procedure  was  executed  under  argon  and  hexane  was 
used  to  wash  the  particles.  The  powder  was  repeatedly  rinsed  with  hexane  until  no  oil  was 
present  on  the  particles. 
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2.1.2  Terfenol-D  composite  manufacture 

The  Terfenol-D  composites  in  this  study  consisted  of  a  polyester  matrix  (Polylite  31610- 
05,  Reichold  Chemicals,  Research  Triangle  Park,  NC)  and  Terfenol-D  particles  prepared 
by  the  procedure  described  in  Section  2.1.1.  The  first  step  of  fabrication  was  to  add  two 
catalysts,  cobalt  naphthenate  (0.1%  by  wt.  resin)  and  N-N  dimenthylaniline  (0.1%  by  wt. 
resin),  and  an  aerator,  Byk-A  501  (0.1%  by  wt.  resin),  to  the  polyester  resin.  After  each  of 
these  was  added,  the  polyester  was  mixed  and  degassed  in  a  vacuum  oven.  Next,  15%  weight 
fraction  (2.2%  volume  fraction)  of  Terfenol-D  particles  was  mixed  into  the  polyester  resin  and 
degassed.  To  initiate  polymerization,  methyl  ethyl  ketone  peroxide  (1.0%  by  wt.  resin)  was 
added.  After  the  final  mixing  and  degassing  step,  the  resin  mixture  was  poured  into  a  silicon 
rubber  mold.  The  mold  was  placed  between  two  permanent  magnets  (Adams  Magnetic 
Products  Co.,  Melrose  Park,  1L)  or  in  an  electromagnet  (Ogallala  Electronics,  Ogallala, 
NE)  and  allowed  to  cure  at  room  temperature  for  48  hours.  The  composite  specimens  were 
12.5  x  12.5  x  25  mm  as  shown  in  Figure  2.1. 

2.1.3  Magnetic  flux  during  cure  and  magnetic  anneal  after  cure 

Applying  a  magnetic  flux  to  a  composite  specimen  during  cure  causes  the  Terfenol-D  particles 
to  orient  into  chains  in  the  direction  of  the  flux.  Each  particle  is  made  up  of  multiple  magnetic 
domains  of  different  orientation.  The  magnetic  flux  causes  the  particles  to  rotate,  such  that 
the  majority  of  the  magnetic  domains  in  a  particle  are  aligned  with  the  magnetic  flux,  and 
agglomerate  into  chains.  When  the  polyester  gels,  the  particles  are  locked  into  this  chained 
configuration.  The  composite  specimens  were  either  cured  with  a  magnetic  flux  through 
their  length  or  a  magnetic  flux  through  their  width  resulting  in  chaining  in  the  direction  of 
the  applied  magnetic  flux  as  shown  in  Figures  2.2  and  2.3. 

Applying  a  magnetic  flux  after  the  specimen  is  fully  cured  will  be  referred  to  in  this 
thesis  as  “magnetic  annealing” .  Magnetic  annealing  partially  orients  the  magnetic  domains 
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in  the  direction  of  the  applied  flux,  but  cannot  rotate  the  entire  particle,  as  the  particles  are 
locked  into  their  position  by  the  solidified  polyester  matrix.  The  specimens  were  magnetically 
annealed  for  20  seconds  either  in  the  electromagnet  or  between  the  two  permanent  magnets 
and  then  tested  immediately  after  the  magnetic  annealing. 

2.2  Experimental  procedure 

All  of  the  specimens  were  tested  in  uniaxial  compression  on  an  Instron  8500  testing  machine 
(Canton,  MA).  LabView™  was  used  to  acquire  the  data  from  the  load  cell,  linear  variable 
displacement  transducer  (crosshead  position),  and  Gaussmeter  (magnetic  flux).  Loading  and 
unloading  was  at  a  constant  rate  of  0.01  kN/s  for  each  test.  The  specimen  was  centered  on 
the  bottom  compression  platen  and  then  an  initial  compressive  load  of  1%  of  the  maximum 
load  was  applied  to  align  the  specimen  before  commencing  the  test.  A  diagram  of  the 
experimental  setup  is  shown  in  Figure  2.4. 

The  compression  platens  were  designed  specifically  for  this  experiment.  They  are  made 
entirely  of  304  austenitic  stainless  steel,  which  is  nonmagnetic.  The  total  length  of  each 
fixture  is  5  inches,  to  separate  the  specimen  from  the  iron  screws  of  the  Instron  machine. 
The  bottom  platen  has  a  spherical  seat  which  helps  to  eliminate  bending  and  off-axis  loading 
when  the  specimen  is  centered  carefully. 

The  magnetic  response  was  measured  with  a  Hall  probe  and  Gaussmeter  (Model  450, 
Lake  Shore  Cryotronics,  Westerville,  OH).  The  Gaussmeter  conditioned  the  signal  from  the 
Hall  probe  and  output  magnetic  flux  (B)  readings.  One  limitation  of  the  setup  was  that 
the  Hall  probe  could  measure  only  one  direction  of  the  magnetic  flux  at  a  time.  To  obtain 
all  three  components  of  the  magnetic  flux  vector,  the  experiment  had  to  be  repeated  three 
times.  Due  to  the  geometry  of  the  probe  as  shown  in  Figure  2.5,  the  magnetic  flux  could 
not  be  measured  at  the  same  distance  from  the  surface  in  all  three  directions.  To  measure 
a  magnetic  flux,  the  probe  had  to  be  oriented  such  that  the  magnetic  flux  went  through  the 
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sensing  area  (a  circle  of  wire).  For  the  measurement,  the  probe  could  be  placed  directly  on 
the  surface  of  the  specimen  and  the  sensing  area  was  approximately  half  the  probe  thickness 
or  0.8  mm  from  the  surface.  For  the  B\  and  B2  measurements,  the  probe  sensing  area  was 
perpendicular  to  the  surface  and  thus  3.8  mm  from  the  surface. 

2.3  Typical  results  and  data  reduction 

For  a  typical  test  of  a  2.2%  Terfenol-D/polyester  composite,  one  component  of  magnetic  flux, 
Bi ,  is  measured  and  plotted  versus  compressive  stress  as  shown  in  Figure  2.6.  The  magnetic 
flux  at  zero  load  will  be  referred  to  as  B™lt.  Since  the  loading  response  is  more  important  than 
the  unloading  response,  the  change  in  magnetic  flux  from  initial  to  maximum  compressive 
stress  is  an  important  quantity  and  will  be  referred  to  as  ABma,x.  In  some  instances,  the 
magnetic  flux  values  were  zeroed  by  subtracting  B™lt.  The  magnetic  flux  is  then  referred  to 
as  the  “zeroed  magnetic  flux”  as  shown  in  Figure  2.7. 

For  the  composite  specimens,  a  special  nomenclature  is  used  to  denote  the  magnetic  flux 
applied  during  cure,  the  annealing  magnetic  flux,  and  the  loading  direction.  Each  test  on  a 
specimen  is  labeled  as  follows:  (a,b,c)  where  a  =  the  direction  of  magnetic  flux  during  cure, 
b  =  the  direction  of  the  annealing  magnetic  flux,  and  c  =  the  loading  direction.  For  example, 
a  test  with  (3,3,1)  indicates  that  both  the  magnetic  flux  during  cure  and  the  annealing 
magnetic  flux  were  in  the  3-direction  while  the  loading  was  in  the  1-clirection.  Further,  the 
graph  axis  is  labeled  as  to  which  component  of  magnetic  flux  ( Bi )  was  measured  in  the 
experiment. 

Since  the  measurements  of  B 1,  Z?2,  and  B3  are  taken  at  different  distances  from  the 
surface  of  the  specimen  as  discussed  in  Section  2.2,  a  study  was  conducted  to  determine  the 
dependence  of  magnetic  flux  measurements  on  the  distance  of  the  probe  sensing  area  from 
the  surface  of  a  specimen.  A  (3,3,1)  Terfenol-D  composite  specimen  was  tested  to  a  uniaxial 
compressive  stress  of  —10  MPa  and  AB giax  was  measured  at  the  midpoint  of  the  specimen  at 
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multiple  distances  up  to  22  mm  from  the  surface  of  the  specimen.  The  annealing  magnetic 
flux  was  —800  G  for  all  tests.  The  results  are  shown  in  Figure  2.8.  It  is  interesting  to  note 
that  while  the  magnetic  flux  dropped  quickly  as  the  probe  was  moved  away  from  the  surface 
of  the  specimen,  the  magnetic  flux  did  not  exhibit  a  1/d3  dependence.  The  magnetic  flux  of 
a  magnetic  dipole  decays  as  1/r3,  so  this  trend  was  expected  [10].  The  dependence  of  the 
magnetic  flux  measurement  on  distance  from  the  specimen  has  a  detrimental  effect  on  the 
results  of  all  three-dimensional  gage  length  characterizations,  since  the  values  of  B i  and  B2 
were  measured  3.8  mm  from  the  surface  while  B%  was  measured  0.8  mm  from  the  surface. 
The  Z?3  measurements  will  be  larger  in  magnitude  than  the  Bi  and  B2  measurements. 

2.4  Gage  length  characterization 

For  the  pure  Terfenol-D  specimen  and  the  Terfenol-D  composites,  measurements  of  magnetic 
flux  were  taken  over  the  gage  length.  The  specimens  were  marked  with  a  grid  as  shown  in 
Figure  2.9.  Then,  magnetic  flux  measurements  were  taken  along  the  middle  of  the  specimen 
at  2  mm  increments  in  the  1-direction  along  the  front  face  perpendicular  to  the  3-direction.  In 
some  cases,  the  side  face  that  is  perpendicular  to  the  —2-direction  was  used  for  measurements. 
Unless  otherwise  noted,  magnetic  flux  was  measured  on  the  front  face. 

A  typical  result  with  the  entire  gage  length  compiled  into  one  graph  is  shown  in  Fig¬ 
ure  2.10.  This  graph  shows  the  variations  in  initial  magnetic  flux  as  well  as  the  shape  and 
magnitude  of  the  loading  and  unloading  responses.  To  compare  the  magnetic  response  in  a 
different  way,  the  magnetic  flux  was  zeroed  by  subtracting  the  magnetic  flux  at  zero  load, 
U]mt,  as  shown  in  Figure  2.11.  This  subtraction  process  removes  the  ambient  magnetic  flux, 
Bfmh ,  of  the  room  as  well  as  the  initial  magnetic  flux  from  the  specimen.  To  measure  the 
ambient  field,  a  polyester  specimen  (nonmagnetic)  identical  in  geometry  and  size  to  the 
composite  specimens  was  loaded  into  the  grips  and  the  three  components  of  the  magnetic 
flux  along  the  gage  length  were  measured.  This  ambient  flux  measurement  procedure  was 
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conducted  before  each  day  of  testing,  and  a  typical  measurement  is  shown  in  Figure  2.12. 
When  the  specimen  was  loaded  there  was  no  change  in  the  magnetic  field.  The  standard 
deviation  from  the  average  magnetic  flux  was  0.001  G.  The  zeroed  response  indicates  how 
much  variation  occurs  in  A5“iax,  but  does  not  reveal  the  direction  of  the  change.  This 
information  is  important  because  it  tells  whether  the  number  of  magnetic  domains  in  the 
direction  of  measurement  is  increasing  or  decreasing. 

2.4.1  Pure  Terfenol-D 

The  pure  Terfenol-D  specimen  was  tested  with  annealing  magnetic  flux  in  the  1-direction 
only.  Annealing  in  any  other  direction  was  not  possible  as  the  specimen  would  rotate  to 
align  its  length  axis  (1-direction)  with  the  magnetic  flux  unless  it  was  physically  restrained 
from  moving.  The  specimen  was  annealed  with  the  permanent  magnets  at  a  magnetic  flux  of 
—800  G  before  each  experiment.  The  specimen  was  loaded  to  —40  MPa  and  then  unloaded. 
The  results  are  included  in  Appendix  A  in  Figures  A.1-A.8.  Measurements  of  B2  were  taken 
on  both  the  front  face  and  the  side  face.  To  compare  the  results,  a  plot  of  A5™ax  versus 
gage  length  position  is  shown  in  Figure  2.13  for  a  uniaxial  compressive  stress  of  —40  MPa. 

In  Figure  2.13,  the  largest  magnetic  flux  change  that  is  constant  across  the  gage  length 
is  A _B™ax.  There  is  a  slight  decrease  in  the  magnetic  flux  measurements  that  were  taken 
closer  to  the  ends  of  the  specimen.  This  decrease  is  also  seen  clearly  in  Figure  A. 2.  Ini¬ 
tially,  the  values  of  B 1  are  between  2  and  5  G,  then  the  flux  decreases  during  loading  to 
approximately  —0.6  G,  which  was  the  ambient  magnetic  flux,  F?amb,  measured  on  that  day 
of  testing.  The  majority  of  the  decrease  in  magnetic  flux  is  in  the  initial  portion  of  the 
loading  (up  to  —5  MPa).  Then,  the  flux  asymptotically  approaches  zero,  indicating  a  satu¬ 
ration  phenomenon.  The  loading  curve  behavior  shows  that  the  domains  have  rotated  away 
from  the  1-direction.  At  —40  MPa,  the  magnetic  flux,  B\  is  nearly  zero  denoting  that  very 
few  domains  remain  in  the  1-direction.  In  the  unloading  portion  of  the  curve,  no  domains 
return  to  their  original  orientation  as  the  B\  values  do  not  change  appreciably.  The  same 
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type  of  saturation  behavior  is  seen  in  the  B2  and  B%  responses.  Most  of  the  magnetic  flux 
change  is  in  the  first  5  MPa  of  loading.  After  that  change,  the  magnetic  flux  asymptotically 
approaches  a  value  with  very  little  change  during  unloading. 

The  changes  in  magnetic  flux,  A 5™ax  and  AB giax,  are  not  constant  across  the  gage 
length  like  AB™*.  Measurements  on  the  front  face  (Figures  A. 3  and  A. 4)  indicate  that 
AB™*  is  very  small  relative  to  the  other  components  of  magnetic  flux,  but  on  the  side  face 
(Figures  A. 5  and  A. 6),  A f?“ax  is  similar  to  A B™x  measured  on  the  front  face  (Figures  A. 7 
and  A. 8).  These  magnetic  response  values  “fan”  across  the  gage  length  with  both  curves 
passing  through  zero  magnetic  flux  near  the  midpoint  of  the  specimen.  The  difference  in  sign 
of  the  two  responses  can  be  explained  by  the  orientation  of  the  probe  during  measurement. 
If  the  response  was  measured  in  a  cylindrical  coordinate  system  where  z  was  along  the  length 
of  the  specimen,  the  2  and  3  axes  would  both  correspond  to  the  r-direction  outward  from  the 
center  of  the  specimen.  In  this  coordinate  system,  the  responses  would  be  nearly  identical 
as  shown  in  Figure  2.14. 

The  “fanning  phenomenon”  can  be  explained  by  examining  the  magnetic  flux  lines  of  a 
uniformly  magnetized  cylinder  shown  in  Figure  2.15.  The  Terfenol-D  specimen  has  a  nearly 
uniform  magnetization  initially.  Most  of  the  domains  are  oriented  in  the  crystal  growth 
direction  and  then  the  annealing  magnetic  flux  orients  even  more  domains  in  the  1-direction. 
The  Terfenol-D  specimen  can  be  modeled  by  the  uniformly  magnetized  cylinder  shown  in 
Figure  2.15.  When  the  magnetic  flux  vector  is  decomposed  into  components  at  two  different 
locations:  one  above  the  midpoint  and  one  below  the  midpoint,  the  “fanning  phenomenon” 
is  observed.  The  Bi  values  are  the  same  above  and  below  the  midpoint,  but  the  B3  values 
are  opposite. 

2.4.2  Terfenol-D  composites 

Terfenol-D  composites  were  created  with  two  different  magnetic  flux  orientations  during  cure. 
The  curing  magnetic  flux  was  either  through  the  thickness  (3-direction)  of  the  specimen,  or 
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through  the  length  (1-direction)  of  the  specimen.  The  magnetic  flux  during  cure  produced 
composites  with  particle  chains  in  the  direction  of  the  applied  flux.  The  strength  of  the  flux 
was  —800  G  for  both  orientations.  Before  each  test,  the  specimen  was  annealed  with  a  field 
of  —800  G  in  the  1,  2,  or  3-direction.  The  composites  were  then  loaded  to  —10  MPa  of 
uniaxial  compressive  stress  and  fully  unloaded. 

2.4.2. 1  Terfenol-D  composites  with  chaining  in  1-direction 

The  gage  length  characterization  results  for  a  (1,1,1)  Terfenol-D  composite  specimen  are 
shown  in  Figures  A.9-A.16.  The  same  Terfenol-D  composite  specimen  was  also  tested  with 
an  annealing  magnetic  flux  in  the  3-direction.  These  (1,3,1)  Terfenol-D  composite  specimen 
results  are  given  in  Figures  A.17-A.24.  Similar  to  the  pure  Terfenol-D  specimen  tests, 
measurements  of  B2  were  taken  on  both  the  front  face  and  the  side  face  for  both  series  of 
tests.  To  compare  the  magnetic  response  at  maximum  load,  Al?“ax,  versus  gage  length 
position  is  shown  in  Figures  2.16  and  2.17  for  a  uniaxial  compressive  stress  of  —10  MPa. 

The  magnetomechanical  behavior  of  the  (1,1,1)  and  (1,3,1)  Terfenol-D  composites  is  sim¬ 
ilar  to  the  pure  Terfenol-D  specimen  discussed  in  Section  2.4.1.  The  value  of  Bi  decreases 
toward  zero  during  loading  for  both  the  composites  and  the  pure  specimen,  but  the  magni¬ 
tude  of  the  magnetic  response  differs.  The  initial  values  of  B i  are  smaller  for  the  composites 
than  the  pure  Terfenol-D  specimen,  indicating  a  difference  in  the  number  of  domains.  While 
the  average  A5“iax  was  —4.1  G  for  the  pure  Terfenol-D  specimen,  A5“iax  averaged  only 
—0.49  G  for  the  (1,1,1)  composite  specimen.  When  the  Terfenol-D  composite  was  magnet¬ 
ically  annealed  in  the  3-direction  ((1,3,1)  composite),  the  average  A5)nax  was  even  smaller, 
—0.17  G.  This  trend  indicates  that  a  larger  Bi  response  is  seen  when  more  domains  are 
oriented  along  the  compressive  stress  direction.  The  pure  Terfenol-D  specimen  had  the  most 
magnetic  domains  that  were  oriented  in  the  1 -direction  due  to  both  manufacturing  and 
annealing.  The  composite  had  only  a  2.2%  volume  fraction  of  Terfenol-D,  so  it  had  a  smaller 
number  of  domains.  When  the  domains  were  oriented  in  the  3-direction  by  annealing,  the 
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response  decreased  even  more  because  some  of  the  domains  were  rotated  away  from  the 
chaining  direction  to  the  direction  of  the  annealing  magnetic  flux. 

On  the  front  face,  AF?2iax  was  nearly  zero  for  both  composites.  The  magnetic  response 
in  the  2  and  3-directions  measured  on  the  side  and  front  faces,  respectively,  showed  the 
“fanning”  phenomenon  seen  in  the  pure  Terfenol-D  specimen.  When  the  magnetic  flux 
changes  are  converted  to  a  cylindrical  coordinate  system,  AB™*  and  AI?“ax  are  closer  in 
magnitude  for  each  composite  as  shown  in  Figure  2.18.  An  interesting  detail  is  that  when 
the  specimen  was  annealed  in  the  3-direction,  the  curves  shifted  up  and  no  longer  passed 
through  zero  at  the  midpoint  of  the  specimen  as  seen  in  the  (1,1,1)  composite  and  pure 
Terfenol-D  specimens. 

The  unloading  response  of  the  composites  does  not  indicate  a  saturation  phenomenon. 
In  Figure  A. 10  corresponding  to  the  (1,1,1)  composite,  the  final  magnetic  flux  (Bi)  after 
unloading  is  about  50%  of  AH™ax.  This  type  of  change  in  magnetic  flux  upon  unloading  was 
seen  in  the  (1,3,1)  composite  as  well. 

2. 4. 2. 2  Terfenol-D  composites  with  chaining  in  3-direction 

A  Terfenol-D  composite  specimen  was  manufactured  with  a  magnetic  flux  applied  through 
its  thickness  (3-direction).  This  specimen  was  then  tested  under  all  three  possible  anneal¬ 
ing  directions.  The  gage  length  characterization  results  for  the  (3,1,1),  (3,2,1)  and  (3,3,1) 
composite  tests  are  given  in  Figures  A.25-A.42. 

The  magnetomechanical  behavior  of  the  composite  chained  in  the  3-direction  varied 
greatly  with  the  direction  of  the  annealing  magnetic  flux.  When  the  specimen  was  annealed 
in  the  1-direction,  as  shown  in  Figure  2.19,  the  magnetic  response  was  quite  similar  to  the 
specimens  discussed  in  Sections  2.4.1  and  2.4.2. 1.  The  average  value  of  AZ?“ax,  —0.078  G,  for 
the  (3,1,1)  specimen  was  lower  than  the  corresponding  values  for  both  the  pure  Terfenol-D 
specimen  and  the  composite  chained  in  the  1-direction.  This  decrease  in  magnitude  is  due 
to  the  initial  alignment  of  domains  in  the  3-direction.  Even  with  the  annealing  magnetic 


21 


flux  in  the  1-direction,  not  all  of  the  domains  are  rotated  to  the  1-direction.  The  final  two 
similarities  are  that  one,  the  fanning  phenomenon  occurred  in  the  B3  response  and  two, 
AB2max  was  approximately  zero. 

For  both  the  (3,2,1)  and  the  (3,3,1)  Terfenol-D  composites,  the  maximum  magnetic  flux 
change  was  in  the  direction  of  the  annealing  field  as  shown  in  Figures  2.20  and  2.21.  The 
value  of  AB 3nax  was  0.42  G  for  the  (3,3,1)  specimen,  while  AB 2liax  was  —0.096  G  for  the 
(3,2,1)  specimen.  The  (3,3,1)  response  was  the  largest  of  all  the  specimens  chained  in  the 
3-direction.  This  behavior  is  similar  to  the  specimens  chained  in  the  1-direction.  The 
responses  are  highest  when  the  direction  of  chaining,  annealing,  and  measurement  is  the 
same.  This  observation  is  further  supported  by  swithcing  the  annealing  magnetic  flnx  from 
the  3-direction  to  the  2-direction  and  seeing  the  average  A F?“ax  drop  an  order  of  magnitude 
from  0.42  G  to  0.040  G. 

The  unloading  responses  of  the  specimen  chained  in  the  3-direction  were  different  for 
each  annealing  direction  as  well.  In  the  (3,1,1)  specimen,  saturation  did  not  occur  and  the 
final  magnetic  flux,  Bi,  was  75%  of  AF?™ax.  Saturation  did  not  appear  to  occur  in  the  (3,2,1) 
or  (3,3,1)  specimens  either,  but  each  had  a  different  unloading  response.  The  final  magnetic 
flux,  B-2  of  the  (3,2,1)  specimen  was  identical  to  Afi™  as  seen  in  Figure  A. 34.  In  contrast, 
the  final  magnetic  flux,  B3 ,  of  the  (3,3,1)  specimen  was  actually  20%  higher  than  A F?™ax. 

2.4.3  Determination  of  piezomagnetic  coefficients 

When  magnetostriction  is  modeled  as  a  first-order  effect,  it  is  called  piezomagnetism,  which  is 
discussed  in  depth  in  Section  1.2.2.  The  governing  equation  for  the  compression  experiments 

is 

Bi  =  dijOj  +  /t ijHj  (2-1) 

where  Bi  denotes  the  magnetic  flnx,  Oj  denotes  the  stresses,  and  <%,-  denotes  the  piezomag¬ 
netic  coefficients.  For  the  compression  experiment,  there  are  no  shear  stresses  and  no  applied 
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magnetic  fields,  which  simplifies  the  governing  equation. 

To  determine  the  piezomagnetic  coefficients  of  the  compression  specimens  studied,  the 
gage  length  results  need  to  be  transformed  to  the  coordinate  system  used  in  the  literature, 
where  the  poling  or  magnetization  axis  is  always  denoted  3.  The  Terfenol-D  specimen  has 
a  magnetization  axis  along  its  length  corresponding  to  its  crystal  growth  axis.  For  the 
Terfenol-D  composites,  the  direction  of  the  curing  held  is  designated  as  the  poling  axis.  With 
these  designations,  there  are  two  different  cases  to  consider.  In  the  first  case,  corresponding 
to  the  pure  Terfenol-D  rod  and  the  (1,*,1)  class  of  composites,  the  poling  axis  (3)  is  along 
the  length  of  the  specimen  and  the  stress  is  also  in  this  3-direction.  In  the  second  case, 
corresponding  to  the  (3,*,1)  composites,  the  poling  axis  is  perpendicular  to  the  loading  axis. 
Both  of  these  geometries  are  shown  in  Figure  2.22. 

In  the  Class  I  experiments,  the  only  nonzero  stress  is  cr3.  The  magnetic  flux  components, 
given  in  the  coordinate  system  used  by  the  literature,  are 


B\  =  di3cr3 

(2.2a) 

B2  =  d23a3 

(2.2b) 

B3  =  d33a3 

(2.2c) 

In  the  Class  II  experiments,  the  only  nonzero  stress  is  Gy. 

given  in  the  coordinate  system  used  by  the  literature,  are 

The  magnetic  flux  components, 

B\  =  d\\(Ji 

(2.3a) 

B2  =  d2\o  1 

(2.3b) 

B3  =  d3i&i 

(2.3c) 

The  numerical  values  for  the  piezomagnetic  coefficients 

were  determined  by  plotting  the 

zeroed  magnetic  flux  versus  stress  for  the  loading  portion  of  an  experiment.  The  response 
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over  the  entire  gage  length  was  used  in  the  linear  curve  fit.  In  the  experiments  where 
the  magnetic  flux  exhibited  the  “fanning”  phenomenon  over  the  gage  length,  the  average 
response  was  nearly  zero,  so  no  piezomagnetic  coefficients  were  calculated. 

2. 4. 3.1  Terfenol-D 

The  magnetic  poling  axis  of  the  Terfenol-D  specimen  is  along  the  length  of  the  specimen  since 
the  specimen  has  a  crystal  growth  axis  along  the  length  and  the  annealing  magnetic  flux  was 
also  along  the  length,  so  the  geometry  is  Class  I  of  Figure  2.22.  The  experimental  data  used 
for  the  curve  fit  are  the  initial  part  (up  to  —2.0  MPa)  of  the  gage  length  characterization  in 
Figure  A. 2.  The  d33  value,  1.48  x  10~10  m/A,  is  lower  by  two  orders  of  magnitude  than  the 
d33  values  reported  for  the  direct  effect  in  Section  1.2.3. 

2. 4. 3. 2  Terfenol-D  composites 

There  were  two  types  of  composites  tested:  those  poled  in  the  direction  of  stress  (Class  I) 
and  those  poled  perpendicular  to  the  stress  (Class  II)  by  curing  with  a  magnetic  flux  strength 
of  800  G.  Within  these  two  classes,  the  direction  of  the  annealing  field  was  varied,  but  its 
magnitude  was  always  800  G.  Only  the  magnetic  responses  that  were  nearly  constant  across 
the  gage  length  were  fit  to  the  piezomagnetic  model.  These  gage  length  loading  responses  for 
the  (1,1,1)  and  (1,3,1)  specimens  are  presented  in  Figures  2.24  and  2.25  with  the  appropriate 
model  predictions.  Similarly,  the  (3,1,1),  (3,2,1),  and  (3,3,1)  specimen  results  are  given  in 
Figures  2.26,  2.27,  and  2.28,  respectively. 

Table  2.1  gives  a  summary  of  the  piezomagnetic  coefficients  determined  for  all  of  the 
Terfenol-D  composites.  The  largest  coefficient  is  d33  =  7.86  x  10~12  m/A  for  the  (1,1,1) 
composite.  When  the  annealing  is  perpendicular  to  the  curing  direction  (the  (1,3,1)  compos¬ 
ite),  the  value  of  d33  decreases  by  almost  80%.  Also,  a  magnetic  flux  in  the  new  annealing 
direction  is  measured  and  yields  another  piezomagnetic  coefficient,  d13,  which  is  even  smaller 
and  negative. 
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Table  2.1:  Comparison  of  piezomagnetic  coefficients  for  the  different  Terfenol-D  composites 
discussed  in  Sections  2.4.2. 1  and  2. 4. 2. 2. 


Type  of 
composite 

Geometry 

class 

Piezomagnetic 
coefficients  (m/A) 

(1.1.1) 

I 

d33  =  7.86  x  10“12 

(1,3,1) 

I 

d33  =  1.73  x  10~12  and  d13  =  -5.48  x  10“13 

(3,1,1) 

II 

du  =  8.37  x  10~ia 

(3,2,1) 

II 

d3l  =  -4.73  x  10"13  and  d2 1  =  1.07  x  10~12 

(3,3,1) 

II 

d31  =  -4.50  x  10~12 

The  largest  coefficient  for  the  Class  II  composites  is  not  as  large  as  0^3  for  the  (1,1,1) 
Tcrfenol-D  composite,  but  it  does  correspond  to  the  specimen  where  the  curing  and  annealing 
directions  are  the  same.  The  largest  Class  II  piezomagnetic  coefficient  is  d3\  with  a  value 
of  —4.50  x  10“12  m/A  for  the  (3,3,1)  composite.  Similar  to  Class  I  phenomenon,  when  the 
annealing  direction  is  perpendicular  to  the  curing  direction,  the  piezomagnetic  coefficient  is 
decreased  in  size  and  new  additional  piezomagnetic  coefficients  are  measured. 

2.5  Effect  of  magnetic  flux  on  cyclic  response 

To  study  the  effects  of  the  magnitude  of  the  curing  and  annealing  magnetic  fluxes  on  the 
Terfenol-D  composites,  specimens  were  manufactured  with  a  curing  flux  of  —2600  G  using  an 
electromagnet.  One  type  of  specimen  was  chained  through  its  thickness  and  denoted  (1,*,1) 
and  one  type  was  cured  through  its  length  and  denoted  (3,*,1).  These  composites  were  not 
tested  for  gage  length  characterization,  but  rather  for  cyclic  response.  Terfenol-D  composites 
created  with  a  curing  magnetic  flux  of  —800  G  were  also  tested  for  cyclic  response. 

To  characterize  the  cyclic  response,  a  specimen  was  annealed  in  its  chaining  direction 
before  testing  in  uniaxial  compression.  For  this  “cycle  1”,  the  specimen  was  loaded  to 
—10  MPa  and  then  completely  unloaded.  For  all  following  cycles,  the  loading  regimen  was 
the  same,  but  the  specimen  was  not  annealed  before  testing.  The  magnetic  flux  was  measured 
in  the  chaining  direction  at  the  midpoint  of  the  front  face  for  each  test.  A  typical  cyclic 
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response  is  shown  in  Figures  2.29  and  2.30.  The  graphs  show  that  the  initial  cycle  is  different 
from  all  following  cycles.  There  is  negligible  difference  between  the  third  cycle  and  the  fiftieth 
cycle,  so  as  few  as  3  cycles  were  necessary  to  observe  a  stable  response. 

Results  for  the  cyclic  response  of  (1,1,1)  Tcrfenol-D  composites  are  shown  in  Figures  B.l- 
B.8.  Similar  results  for  (3,3,1)  Terfenol-D  composites  are  given  in  Figures  B.9-B.14.  The 
curing  flux  and  annealing  flux  were  either  —800  G  or  —2600  G  for  each  test.  One  interesting 
trend  in  all  the  plots  of  nonzeroed  magnetic  flux  is  that  the  unloading  response  of  cycle  1  is 
the  same  as  the  loading  response  of  the  following  cycle.  Then,  the  unloading  of  that  cycle 
traces  back  the  same  unloading  response  as  cycle  1.  This  is  repeated  in  all  later  cycles,  there 
is  little  difference  between  any  cycles  after  the  first.  The  trend  shows  that  the  annealing 
held  aligns  some  domains  in  the  direction  of  chaining,  but  these  domains  do  not  return  to 
the  chaining  direction  upon  unloading.  After  the  first  cycle,  the  domains  that  were  locked 
into  position  by  the  curing  flux  are  solely  responsible  for  the  measured  magnetic  response. 
Even  without  any  annealing,  a  repeatable  measurable  response  is  seen.  The  nonannealed 
response  is  smaller  than  the  first  cycle  response  with  an  annealing  magnetic  flux. 

By  varying  the  magnitude  of  the  curing  and  annealing  magnetic  fluxes,  the  initial  cycle 
response  and  final  cycle  response  can  be  changed.  Table  2.2  shows  a  comparison  of  AR™ax  for 
(1,1,1)  Terfenol-D  composites.  Increasing  the  curing  magnetic  flux  from  —800  G  to  —2600  G 
increased  the  values  of  AR™ax  for  both  the  first  cycle  and  the  final  cycle.  Increasing  the 
annealing  magnetic  flux  from  —800  G  to  —2600  G  increased  the  initial  cycle  response,  but 
did  not  have  a  significant  effect  on  the  final  cycle  response. 

The  dependence  of  A R™ax  on  curing  and  annealing  magnetic  flux  for  (3,3,1)  Tcrfenol-D 
composites  presented  in  Table  2.3.  In  these  composites,  increasing  the  curing  flux  actually 
lowered  the  magnitude  of  AB in  the  first  cycle.  Conversely,  increasing  the  curing  flux 
increased  the  final  cycle  response.  As  seen  in  the  (1,1,1)  composites,  the  annealing  flux 
strength  had  very  little  effect  on  the  final  cycle  response  of  the  (3,3,1)  composites. 
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Table  2.2:  Comparison  of  Af?“iax  for  different  combinations  of  curing  and  annealing  magnetic 
fluxes  for  (1,1,1)  Terfenol-D  composites. 


Curing 

magnetic  flux  (G) 

Annealing 
magnetic  flux  (G) 

ASf8*  (G) 

First  cycle 

Final  cycle 

-800 

-800 

0.55 

0.30 

-800 

-2600 

0.68 

0.34 

-2600 

-800 

0.75 

0.45 

-2600 

-2600 

0.77 

0.42 

Table  2.3:  Comparison  of  AB1™*  for  different  combinations  of  curing  and  annealing  magnetic 
fluxes  for  (3,3,1)  Terfenol-D  composites. 


Curing 

magnetic  flux  (G) 

Annealing 
magnetic  flux  (G) 

A B“ax  (G) 

First  cycle 

Final  cycle 

-800 

-800 

0.44 

-0.04 

-2600 

-800 

-0.16 

-0.82 

-2600 

-2600 

-0.20 

-0.89 

Figure  2.1:  Geometry  of  uniaxial  compression  specimen. 
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Magnetic  flux 
applied  during  cure 


-<  0.92  mm 


Entire  front  face 


Magnified  image 


Figure  2.2:  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  with 
a  magnetic  flux  through  the  length. 


Entire  front  face 


Magnified  image 


Figure  2.3:  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  with 
a  magnetic  flux  through  the  width. 


Figure  2.4:  Experimental  setup  for  uniaxial  compression  tests. 


Figure  2.5:  Schematic  of  Hall  probe  used  to  measure  magnetic  flux. 
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Magnetic  flux,  B  (G) 


0  -2  -4  -6  -8 

Compressive  stress,  cM  (MPa) 

Figure  2.6:  Typical  response  of  (1,1,1)  composite  specimen  under  uniaxial  com 
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Compressive  stress,  cn  (MPa) 


Figure  2.7:  Typical  response  of  (1,1,1)  composite  specimen  under  uniaxial  compressive  stress 
with  magnetic  flux  zeroed  by  subtracting  magnetic  flux  value  at  zero  load, 
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Magnetic  flux  change  over  -10  MPa,  AB"lilx  (G) 


Figure  2.8:  Distance  dependence  of  AD™ax  measurement  on  a  (3,3,1)  Terfenol-D  composite 
specimen  loaded  to  —10  MPa  and  annealed  at  —800  G  before  each  test. 


25  mm 


Figure  2.9:  Compression  specimen  marked  with  grid  for  gage  length  characterization. 
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Compressive  stress,  on  (MPa) 


Figure  2.11:  Typical  gage  length  response  of  a  (1,1,1)  Terfenol-D  composite  with  magnetic 
flux  zeroed  by  subtracting  B\mt. 


Figure  2.12:  Typical  ambient  magnetic  flux  measured  on  a  polyester  specimen. 
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Figure  2.13:  Comparison  of  ABmax ,  the  change  in  magnetic  flux  over  —40  MPa,  for  a  pure 
Terfenol-D  specimen  annealed  at  a  magnetic  flux  of  —800  G  in  the  1-direction. 
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Position  along  gage  length  (mm) 


Figure  2.16:  Comparison  of  AZ?max,  the  change  in  magnetic  flux  over  —10  MPa,  for  a  (1,1,1) 
Terfenol-D  composite  specimen  cured  and  annealed  at  —800  G. 
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Magnetic  flux  change  over  -10  MPa,  ABmax  (G) 


Magnetic  flux  change  over  -10  MPa,  ABmax  (G) 


Position  along  gage  length  (mm) 


Figure  2.18:  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  A5max,  for  (1,1,1) 
and  (1,3,1)  Terfenol-D  composite  specimens  cured  and  annealed  at  a  magnetic  flux  of  —800  G, 
where  z  is  the  axis  along  the  length  of  the  specimen  and  r  points  outward. 
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Figure  2.19:  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  a  (3,1,1) 
Tcrfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G. 


Magnetic  flux  change  over  -10  MPa,  ABmax  (G) 
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Figure  2.20:  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  ABmax,  for  a  (3,2,1) 
Tcrfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G. 


Position  along  gage  length  (mm) 


Figure  2.21:  Comparison  of  the  change  in  magnetic  flux  over  —10  MPa,  AHmax,  for  a  (3,3,1) 
Tcrfenol-D  composite  specimen  cured  at  —800  G  and  annealed  at  —800  G. 


Class  I:  Poling  axis  is  the  same  Class  II:  Poling  axis  is  perpendicular 

as  the  loading  axis  to  the  loading  axis 

Figure  2.22:  Geometries  of  compression  specimens  converted  to  coordinate  system  used  in 
piezomagnetic  literature. 
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Figure  2.23:  Piezomagnetic  model  prediction  and  experimental  data  for  Terfenol-D  specimen 
poled  in  the  same  direction  as  the  compressive  stress. 
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Zeroed  magnetic  flux  in  poling  direction,  B,  (G) 


Compressive  stress,  o3  (MPa) 


Figure  2.24:  Piezomagnetic  model  prediction  and  experimental  data  for  (1,1,1)  Tcrfenol-D 
composite  specimen. 


Zeroed  magnetic  flux,  B  (G) 
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Figure  2.25:  Piezomagnetic  model  prediction  and  experimental  data  for  (1,3,1)  Tcrfenol-D 
composite  specimen. 
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Compressive  stress,  Oj  (MPa) 


Figure  2.26:  Piezomagnetic  model  prediction  and  experimental  data  for  (3,1,1)  Tcrfenol-D 
composite  specimen. 
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Figure  2.27:  Piezomagnetic  model  prediction  and  experimental  data  for  (3,2,1)  Tcrfenol-D 
composite  specimen. 
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Compressive  stress,  Oj  (MPa) 

Figure  2.28:  Piezomagnetic  model  prediction  and  experimental  data  for  (3,3,1)  Terfenol-D 
composite  specimen. 
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Chapter  3 

Residual  stress  analysis 


3.1  Introduction 

Several  competing  mechanisms  contribute  to  residual  stress  development  during  the  fab¬ 
rication  of  Terfenol-D  composites.  Initially,  the  polyester-particulate  mixture  is  at  room 
temperature  and  able  to  flow  freely.  Then,  the  mixture  is  poured  into  a  mold  and  placed 
between  two  magnets.  The  magnetic  held  causes  the  Terfenol-D  particles  to  strain  via  the 
magnetostrictive  effect  as  discussed  in  Section  1.2.  At  this  time,  the  polyester  is  also  gener¬ 
ating  heat  from  the  polymerization  and  crosslinking  reactions.  Until  the  polymer  “sets”  or 
“gels” ,  there  is  no  stress  on  the  particles  and  they  move  freely  to  align  in  chains  along  the 
held  lines  between  the  magnets  as  shown  in  Figure  3.1.  Once  gelation  occurs,  the  particles 
are  constrained  and  any  further  chemical  shrinkage  of  the  matrix  will  produce  stresses  on  the 
particles.  Also,  the  temperature  drops  slowly  back  down  to  room  temperature  and  thermal 
mismatch  stresses  occur.  When  a  specimen  is  taken  out  of  its  mold,  it  is  no  longer  under 
the  influence  of  a  magnetic  held,  so  the  magnetostrictive  strains  are  released. 

Many  theories  and  models  exist  for  predicting  residual  stresses  in  a  composite,  but  few,  if 
any,  can  be  directly  applied  to  the  particulate  composites  at  hand.  Some  of  the  papers  focus 
on  hber-reinforced  composites,  especially  with  thermosetting  resin  matrices  [34,  31].  These 
theories  could  be  applied  to  the  Terfenol-D  particulate  composites  at  hand  if  the  particulate 
chains  are  modeled  as  hbers.  There  are  also  a  large  number  of  papers  in  the  literature  for 
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metal-matrix  particulate  composites  [3,  12,  24],  The  major  focus  of  these  papers  is  residual 
stresses  induced  by  thermal  mismatch.  These  papers  offer  the  best  theories  to  build  on  and 
modify  for  Terfenol-D  composites. 

3.2  Analysis 

A  modified  Eshelby  model  is  used  to  determine  the  residual  stress  field  in  the  Terfenol-D 
composite.  The  model  is  based  on  a  similar  analysis  presented  by  Hsueh  and  Becher  [12] 
in  which  ellipsoidal  inclusions  are  modeled.  The  residual  stress  from  thermal  mismatch  is 
calculated  for  a  number  of  different  inclusion  geometries  using  an  anisotropic  coefficient  of 
thermal  expansion.  The  geometry  of  the  inclusion  is  represented  by 


with  A3  being  the  unique  axis  of  the  ellipsoid  as  shown  in  Figure  3.2.  By  choosing  the  ratio 
c/a  appropriately,  the  geometry  of  an  ellipsoid,  sphere,  disk,  or  fiber  can  be  represented.  For 
this  analysis,  the  inclusions  are  modeled  in  two  ways.  First,  a  single  particle  is  modeled  as 
a  spherical  inclusion.  This  may  not  be  the  best  geometric  assumption  as  the  particles  group 
very  close  together  and  form  chains  as  shown  in  Figure  3.1.  Therefore,  in  the  second  analysis, 
a  chain  of  multiple  particles  is  represented  by  an  ellipsoidal  inclusion.  Linear  isotropic 
elasticity  is  assumed  for  both  the  matrix  and  the  inclusions.  The  model  is  expanded  from 
the  work  of  Hsueh  and  Becher  [12]  to  include  magnetostrictive  effects  and  chemical  shrinkage 
as  sources  of  residual  stress. 

The  Eshelby  model  represents  the  problem  of  mismatch  stresses  by  applying  a  series 
of  virtual  operations  on  an  inclusion.  First,  it  is  imagined  that  the  inclusion  is  removed 
from  the  matrix  and  all  mismatch  strains  are  allowed  to  develop  unconstrained.  The  strains 
that  develop  in  this  step  are  defined  as  the  transformation  strains,  e **.  In  this  analysis,  the 
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transformation  strains  are 


e%  =  +  4  +  e?4  (3.2) 

where  eff  are  the  thermal  mismatch  strains,  eff  are  the  chemical  shrinkage  strains,  and  e™ag 
are  the  magnetostrictive  strains.  All  these  strains  are  either  isotropic  or  transversely  isotropic 
about  the  X3  axis  of  the  inclusion. 

Next,  an  imaginary  surface  traction  is  applied  to  the  inclusion  to  fit  it  back  into  its  void 
in  the  matrix.  The  surface  traction  creates  a  stress  in  the  inclusion  but  no  stress  in  the 
matrix.  To  correct  this  imbalance,  the  applied  surface  traction  is  treated  as  a  layer  of  body 
force  and  an  equal  and  opposite  layer  of  body  force  is  applied  at  the  interface  to  make  the 
net  surface  traction  vanish  along  the  interface.  The  strains  induced  by  these  body  forces  are 
termed  constrained  strains,  ef;,  and  are  functions  of  the  transformation  strains,  e**.  The  final 
residual  stresses  in  the  inclusion  are  spatially  uniform  and  arise  from  the  difference,  eci3  —  e(* , 
while  the  residual  stresses  in  the  matrix  arise  from  ef,-. 

The  Eshelby  model,  as  presented  for  an  ellipsoidal  inclusion  in  an  infinite  matrix,  gives 
the  relationship  between  e-  and  ej*  for  the  case  where  the  clastic  properties  of  the  matrix 
and  inclusion  are  the  same.  For  a  composite  consisting  of  two  distinct  phases,  an  additional 
modification  is  required.  First,  isotropic  linear  elasticity  is  assumed  for  both  the  matrix 
and  the  inclusion.  The  true  inclusion  is  then  virtually  replaced  by  an  equivalent  inclusion 
and  the  equivalent  transformation  strains,  eh,  result.  To  insure  the  sameness  of  stresses  and 
displacements  in  the  two  inclusions,  two  relationships  must  hold: 

A'e(ec  —  e**)  =  Kjf-J)  (3.3a) 

—  £11)  =  —  e^)  (3.3b) 

where  K  and  G  are  bulk  and  shear  moduli  for  the  ellipsoidal  inclusion  (subscript  e)  and  the 
matrix  (subscript  m).  The  dilatational  and  deviatoric  strains  in  Equations  (3.3a)  and  (3.3b) 
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are  defined  by 


e  —  en  +  e 22  +  633  (3.4a) 

1 

en  =  -  (en  —  £33)  (3.4b) 

Finally,  the  Eshelby  model  is  modified  to  reflect  the  finite  volume  fraction  of  inclusions. 
The  presence  of  surrounding  inclusions  introduces  an  additional  constraint,  such  that 


e 


C 

V 


6  ■  +  6  • 
V  1  V 


(3.5) 


where  ec'-  are  the  constrained  strains  in  an  infinite  matrix  and  ef"  are  the  additional  strains 

LJ  LJ 

of  the  surrounding  inclusions.  The  first  constrained  strains  are 


e ij  ^ijkl^kl 


(3.6) 


where  Sijki  is  the  Eshelby  tensor.  The  second  constrained  strains  are  calculated  from  me¬ 
chanical  equilibrium  to  give 


“u 


(3.7) 


where  /  represents  the  volume  fraction  of  inclusions.  Finally,  the  total  constrained  strains 
are 

eij  =  (1  -  /)  SijklAl  +  felj  (3-8) 
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The  Eshelby  tensor  coefficients,  Sl3kh  for  ellipsoidal  inclusions  are 


S3333  —  Qc2  ICc  +  RIc 

'S'llll  =  *5*2222  =  Qd2Iaa  +  RI a 

*5*1133  =  *S*2233  =  Qc2 1 ca.  ~  RI a 

*S*1122  =  *S*2211  =  Qu2 1  ah  ~~  RI a 

*S*3311  =  *S*3322  =  Q(2Ica,  —  RIc 


where  Q  and  R  are  given  by 


8vr  (1  -  vm) 


8vr  (1  -  vm) 


(3.9a) 

(3.9b) 

(3.9c) 

(3.9d) 

(3.9e) 


(3.10a) 

(3.10b) 


When  the  particle  chains  are  modeled  as  an  ellipsoid,  the  coefficients  are  given  by 


L  = 


2ira2c 


(c2  —  a2)2 


-  Wi 

a  \  or 


Ir  =  47T  —  2/a 


cosh  1  (  - 


Ira  Icb 


h-Ic 

3  (c2  -  a2 


I rr 


lab  — 


I aa 


--2I 

3  c2 
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3  a2 
3/ab 


-‘ca 

T 


(3.11a) 

(3.11b) 

(3.11c) 

(3. lid) 
(3. lie) 
(3. Ilf) 


For  the  spherical  inclusions,  the  coefficients  are  given  by 


47T 

T 


3/cb 


3/ab 


47T 
5 a2 


(3.12a) 

(3.12b) 
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The  equivalent  transformation  strains  are  found  by  solving  Equations  (3.3a),  (3.3b),  and 


(3.8).  For  a  spherical  inclusion, 


-li 


=  e 


22 


“33 


(2 AKe  +  BGe)  eft  +  (AKe  -  BGe)  ef3 
3  AB 

(2 AKe  -  2 BGe)  e\\  +  ( AKe  +  2 BGe)  e|*3 
3AB 


(3.13a) 

(3.13b) 


where  the  constants  A  and  B  are  as  follows: 


A  —  Gm  +  (/  +  (1  —  /)  (S3333  —  S1133))  (Ge  —  Gm)  (3.14a) 

B  =  Am  +  (/ +  (1  — /)  (53333  +  21S1133))  {Ke  —  Am)  (3.14b) 


The  equivalent  transformation  strains  for  an  ellipsoidal  inclusion  with  X3  as  the  major  axis 
are 


Ai 


“33 


(2 AKe  +  BGe )  e(j  +  (AKe  -  BGe )  e& 
CA-BD 

{2DKe  +  CGe )  e‘j  +  (DKe  -  CGe )  e% 
BD  —  CA 


(3.15a) 

(3.15b) 


where  the  constants  A  and  B  are  defined  in  Equations  (3.14a)  and  (3.14b),  and  C  and  D 
are  as  follows: 


C  =  2Am  +  2  (/  +  (1  -  /)  (53311  +  Sun  +  Si  122))  {Ke  -  Km)  (3.16a) 

D  =  —Gm  +  (— /  +  (1  —  /)  (2S3311  —  Aim  —  51122))  (Ge  —  Gm)  (3.16b) 


The  constrained  strains,  efj,  are  calculated  from  the  equivalent  transformation  strains  via 
Equation  (3.8).  Finally,  using  the  stress-strain  relation,  the  residual  stresses  in  an  inclusion 
are 


O'ij  — 


v  E 


(1  +  z/m)  (1  —  2um) 


(ec  -  el)  Sij  + 


E„ 


1  +  vn 


(ec-  -  e--) 
\  >3  y ) 


(3.17) 
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Table  3.1:  Summary  of  material  properties  used  in  modified  Eshelby  model  for  residual 
stresses. 


E  [GPa] 

V 

a  [“C"1] 

/ 

Polyester 

3.8 

0.33 

75  x  lO”6 

0.9776 

Terfenol-D 

50 

0.444 

12  x  10~6 

0.0224 

3.3  Results 

3.3.1  Material  properties 

In  this  section,  the  modified  Eshelby  model  of  Section  3.2  is  used  to  analyze  residual  stresses 
in  the  inclusions  in  the  composites  described  in  Section  2.1.2.  The  necessary  material  prop¬ 
erties  and  other  parameters  for  the  Terfenol-D  particulate  composites  were  found  in  the 
literature  [32]  and  from  the  suppliers.  A  summary  of  the  material  properties  is  presented  in 
Table  3.1.  The  Terfenol-D  particles  are  20-40  //rn  in  their  longest  dimension,  so  a  diameter 
of  30  /ini  is  assumed  when  modeling  a  single  particle  as  a  spherical  inclusion.  A  typical  com¬ 
posite  specimen  was  studied  under  an  optical  microscope  and  fifty  chains  of  particles  were 
measured.  The  average  chain  length  in  the  3-direction  was  955  ±  469  /irn  and  the  average 
width  was  110  ±  33  /irn.  In  the  ellipsoidal  particle  chain  model,  the  constants  a  and  c  are 
then  55  /rm  and  478  /mi,  respectively.  The  model  geometries  are  shown  in  Figure  3.3. 

3.3.2  Thermal  mismatch 

The  Terfenol-D  composites  in  this  study  have  a  polyester  matrix  that  is  cured  at  room 
temperature.  As  the  polyester  cures,  heat  is  released  in  the  polymerization  reaction.  Initially, 
the  polymer  is  a  liquid  with  low  viscosity,  and  the  particles  in  the  composite  are  assumed  to  be 
stress-free.  Then,  at  a  point  known  as  the  gel  point,  the  polymer  sets  and  begins  to  solidify. 
When  gelation  occurs  the  polymer  is  at  an  elevated  temperature  due  to  the  heat  created 
by  the  polymerization  reaction.  The  temperature  at  this  point  is  known  as  the  stress-free 
temperature.  As  the  composite  cools,  thermal  residual  stresses  develop.  These  stresses  arise 
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from  the  difference  in  the  coefficients  of  thermal  expansion  of  the  two  constituent  phases.  If 
the  assumption  of  isotropic  thermal  expansion  is  made,  the  transformation  strains  for  the 
Eshelby  model  presented  in  Section  3.2  are 

4*1  =  4*2  =  4  =  («e  -  am)  AT  (3.18) 

where  AT  is  the  final  temperature  minus  the  stress-free  temperature.  The  residual  stresses  in 
the  inclusions  are  compressive  and  increase  with  increasing  magnitude  of  AT.  The  residual 
stresses  on  a  spherical  particle  are  hydrostatic,  while  the  residual  stresses  on  an  ellipsoidal 
particle  chain  are  biased  in  the  chain  direction  as  shown  in  Figure  3.4. 

To  determine  the  stress-free  temperature  and  the  gel  point  of  the  Terfenol-D  composites, 
a  simple  experiment  was  performed.  A  polyester  specimen  was  manufactured  and  a  type  K 
thermocouple  was  embedded  in  the  center  of  the  specimen  during  the  cure  cycle.  All  other 
steps  of  the  manufacturing  process  were  exactly  the  same  as  those  described  in  Section  2.1.2 
except  that  no  particles  were  mixed  into  the  resin.  Lab  VIEW™  was  used  to  acquire  the 
temperature  data.  Multiple  measurements  were  taken,  averaged,  and  then  output  to  a  hie 
every  15  seconds.  An  identical  specimen  was  manufactured  simultaneously,  but  not  placed 
between  the  magnets  or  equipped  with  a  thermocouple.  The  gel  point  was  recorded  as  the 
time  when  a  probe  could  no  longer  penetrate  the  surface  of  this  specimen.  The  results  from 
the  test  showed  the  gel  point  to  be  50  minutes  corresponding  to  a  stress-free  temperature  of 
48°C  as  shown  in  Figure  3.5.  The  change  in  temperature,  AT,  in  Equation  (3.18)  is 

AT  =  24° C  -  48°C  =  -24°C 

The  thermally  induced  residual  stresses  from  this  AT  for  a  single  spherical  inclusion  are 

oh1)  =  (Too  =  (Too  =  —8.0  MPa 
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and  the  residual  stress  for  a  particle  chain  are 


_th 


_  th 

—  (J22 


<7 


th 

33 


-5.7  MPa 
-38.8  MPa 


3.3.3  Chemical  shrinkage 

Another  effect  of  the  polymerization  reaction  is  chemical  shrinkage.  As  the  monomers  join 
to  form  polymer  chains  and  crosslinking  occurs  between  the  chains,  the  volume  of  polymer 
decreases.  The  chemical  shrinkage  can  occur  both  before  and  after  gelation  and  can  be  as 
large  as  20%  [28].  Before  gelation,  the  shrinkage  causes  no  stress  on  the  Terfenol-D  particles 
because  the  polymer  relaxes  immediately,  but  after  gelation  the  shrinkage  compresses  the 
particles  causing  the  development  of  residual  stresses.  The  strains  associated  with  volumetric 
chemical  shrinkage  are  isotropic  and  calculated  via 

Dl  =  ^22  =  ^33  =  (1  +  AK.h)3  —  1  (3.19) 

where  AV%  is  the  volumetric  chemical  shrinkage  [34],  The  transformation  strains  for  the 
Eshelby  model  are  the  shrinkage  strains  of  the  inclusion  minus  the  shrinkage  strains  of  the 
matrix.  Since  there  is  no  chemical  shrinkage  of  the  inclusion,  the  transformation  strains  are 

4  =  4*2  =  4  =  -  ((1  +  AKh)»  -  l)  (3.20) 

The  residual  stress  state  clue  to  chemical  shrinkage  is  hydrostatic  pressure  on  the  Terfenol-D 
particles  modeled  as  spherical  inclusions.  The  pressure  increases  with  higher  amounts  of 
shrinkage  as  shown  in  Figure  3.6.  For  the  ellipsoidal  particle  chain,  the  residual  stresses  are 
negative  as  well,  but  0*33  is  almost  an  order  of  magnitude  larger  than  a\\. 

For  the  Terfenol-D  composites  in  this  study,  an  experiment  was  conducted  to  determine 
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Table  3.2:  Summary  of  chemical  shrinkage  experiment  results. 


Specimen  # 

«  [g] 

b  [g] 

c  [g] 

SG 

Vf  [cm3] 

Vi[cm3} 

AVch 

1 

4.8882 

0.4649 

1.6946 

1.3361 

3.667 

3.906 

-6.12% 

2 

4.9054 

0.4649 

1.7053 

1.3384 

3.674 

3.906 

-5.95% 

the  total  volumetric  chemical  shrinkage.  An  upper  bound  to  the  shrinkage  was  determined  in 
the  experiment,  but  this  bound  included  both  the  shrinkage  before  gelation  and  the  shrinkage 
afterward.  The  experiment  was  conducted  according  to  ASTM  standard  D  792-91.  First,  the 
specimen  was  weighed  in  air  and  its  weight  (a)  was  recorded.  Then,  a  wire  was  suspended 
from  the  scale,  immersed  in  water,  and  weighed  (6).  Finally,  the  specimen  was  suspended 
by  the  wire,  immersed  in  the  water,  and  its  weight  was  recorded  (c).  The  specific  gravity 
was  calculated  via 


SG 


a 

a  +  c  —  b 


(3.21) 


and  the  final  volume,  Vf,  was 


Vi 


a 


P  water  S  G 


(3.22) 


where  pwater  was  the  density  of  water  at  room  temperature  (997.6  kg/m3  at  22.5°C).  The  ini¬ 
tial  volume,  Vi,  was  determined  by  measuring  the  size  of  the  specimen  mold.  The  volumetric 
chemical  shrinkage  is 

Vf  —  V 

AVch  =  (3.23) 


Two  specimens  were  tested  and  the  results  are  shown  in  Table  3.2.  The  average  volumetric 
chemical  shrinkage  was  —6.0%.  Using  the  estimate  that  only  half  of  this  shrinkage  (—3.0%) 
occurs  after  gelation,  the  chemically  induced  residual  stresses  on  a  single  spherical  particle 
are 

df!  =  an,  =  (Too  =  —53.5  MPa 
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and  for  a  particle  chain  are 


ch 

^11 


_  _ch 

—  a22 


a ■ 


ch 

33 


-37.8  MPa 
-258.9  MPa 


3.3.4  Magnetostrictive  effects 

The  final  mechanism  that  creates  residual  stresses  is  the  magnetostrictive  effect.  Since  the 
specimens  are  cured  between  two  permanent  magnets  as  discussed  in  Section  2.1.2,  the 
Terfenol-D  particles  strain  in  response  to  the  magnetic  field,  up  to  a  possible  saturation 
magnetostriction  strain  of  2.4  x  10~3  in  the  direction  of  the  held.  An  assumption  of 


mag  _  mag  _  1  mag 

611  ~  e22  ~  ~~  2  633 


(3.24) 


is  made  for  the  strains  in  the  plane  of  isotropy.  These  strains  are  then  used  as  the  trans¬ 
formation  strains  in  the  Eshelby  model,  but  must  be  negative  as  they  are  released  after  the 
composite  is  removed  from  the  magnetic  held: 


J*  _  _  mag 

fc33  —  fc33 

1 

t*  _  t*  _  _  mag 

ell  —  e22  —  2 633 


(3.25a) 

(3.25b) 


The  magnetostrictive  effect  on  residual  stresses  is  shown  in  Figure  3.7.  In  the  direction  of 
the  applied  held  (3),  the  residual  stress  is  positive,  but  in  the  plane  of  isotropy,  the  residual 
stresses  are  negative  and  smaller.  For  the  manufacturing  process  described  in  Section  2.1.2, 
the  magnetic  held  on  the  composites  is  900  Oe,  which  corresponds  to  e^ag  =  0.55  x  10~3 
as  shown  in  Figure  3.8  [21]  for  a  Terfenol-D  rod  with  no  prestress.  These  values  of  magne- 
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tostrictive  strains  correspond  to  residual  stresses  of 


mag  mag 

°n  =  ^22 


<7 


mag 

33 


-0.8  MPa 
1.6  MPa 


for  the  spherical  particle  model  and 


a™s  =  a™e 


a- 


mag 

33 


-0.2  MPa 
12.5  MPa 


for  the  ellipsoidal  particle  chain  model  when  only  magnetostrictive  effects  are  considered. 


3.3.5  Combination  of  all  effects 

Combining  all  of  the  residual  stress  mechanisms  discussed  in  Sections  3.3.2,  3.3.3,  and  3.3.4 
gives  transformation  strains  of 

=  (ae  —  am)  AT  —  (1  +  AVdi)®  +  1  +  — e33  8  (3.26a) 

=  (ae  -  am)  AT  -  (1  +  AVch) ^  +  1  -  e“3ag  (3.26b) 

for  the  Eshelby  model  discussed  in  Section  3.2.  Using  the  results  from  the  experiments  on  the 
Terfenol-D  composites  in  Sections  3.3.2,  3.3.3,  and  3.3.4,  the  total  residual  stresses  and  the 
contribution  from  each  mechanism  are  summarized  in  Table  3.3.  Chemical  shrinkage  has  the 
largest  effect  on  the  magnitude  of  the  residual  stresses  for  both  models.  The  magnetostrictive 
effect  gives  the  smallest  contribution,  and  is  distinct  in  that  <733  due  to  the  magnetostrictive 
effect  is  positive  whereas  all  other  residual  stress  contributions  are  negative. 


ct* 

-11 


=  e 


22 


ct* 

“33 
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Table  3.3:  Summary  of  residual  stresses  for  spherical  particle  model  and  ellipsoidal  particle 
chain  model. 


spherical  particle  model 

particle  chain  model 

<Tii  (MPa) 

0-33  (MPa) 

(T11  (MPa) 

0-33  (MPa) 

Thermal  mismatch  (AT  =  — 24°C) 

-8.0 

-8.0 

-5.7 

-38.8 

Chemical  shrinkage  (AVA  =  —3.0%) 

-53.5 

-53.5 

-37.8 

-258.9 

Magnetostriction  (e^ag  =  0.00055) 

-0.8 

1.6 

-0.2 

12.5 

Total  residual  stress 

-62.3 

-59.9 

-43.7 

-285.2 

3.4  Summary 


The  residual  stresses  in  Terfenol-D  particulate  composites  are  derived  using  a  modified 
Eshelby  model.  Closed-form  analytical  solutions  are  obtained  for  a  spherical  particle  and 
a  chain  of  particles  modeled  as  an  ellipsoid.  Linear  isotropic  elasticity  is  assumed  for  the 
matrix  and  the  inclusions.  The  residual  stresses  are  caused  by  three  mechanisms:  thermal 
mismatch,  chemical  shrinkage  of  the  matrix,  and  magnetostriction  of  the  particles.  Thermal 
mismatch  and  chemical  shrinkage  cause  negative  residual  stresses  (Figures  3.4  and  3.6),  while 
the  magnetostrictive  strains  cause  a  positive  residual  stress  in  the  direction  of  the  applied 
held  (3)  and  negative  stresses  in  the  plane  of  isotropy  perpendicular  to  the  held  (Figure  3.7). 

Several  experiments  were  conducted  to  determine  the  important  parameters  for  the 
Terfenol-D  particulate  composites  in  this  study.  The  stress-free  temperature  of  the  com¬ 
posite  was  measured  to  be  48°C  (Figure  3.5).  The  total  chemical  shrinkage  was  measured  to 
be  —6.0%,  and  half  of  that  shrinkage  (—3.0%)  was  assumed  to  occur  after  gelation.  The  mag¬ 
netic  held  during  cure  was  measured  to  be  900  Oe,  leading  to  e^ag  =  0.00055  (Figure  3.8).  If 
the  assumptions  that  the  magnetostrictive  strains  in  the  plane  of  isotropy  perpendicular  to 
the  applied  held  are  half  of  e^ag  and  the  volume  fraction  of  Terfenol-D  is  0.0224,  the  residual 
stresses  are 


=  22  =  —62.3  MPa 

(Too*  =  —59.9  MPa 
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for  a  single  particle  modeled  as  a  spherical  inclusion  and 


=  C22  =  -43.9  MPa 
atot  =  _286.8  MPa 


for  a  chain  of  multiple  particles  modeled  as  an  ellipsoidal  inclusion. 


Xj 


2.75  mm - ► 


n 


H  (magnetic  field 
applied  during  cure) 


Figure  3.1:  Particle  chaining  phenomenon  in  Terfenol-D  particulate  composites  cured  under 
a  magnetic  field  in  X3  direction  [2], 
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X, 


X, 


Figure  3.2:  Geometry  of  the  ellipsoidal  inclusion  in  the  modified  Eshelby  model. 
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Particle  chain  modeled  as  a 
single  ellipsoidal  inclusion 


Figure  3.3:  Model  geometries  for  the  modified  Eshelby  model. 
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Final  temperature  minus  stress-free  temperature,  AT  (  C) 


Figure  3.4:  Residual  stresses  from  thermal  mismatch. 


Figure  3.6:  Residual  stresses  from  chemical  shrinkage 


Figure  3.7:  Residual  stresses 


0.0016 


0.0020 
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magnetostrictive  strains. 


Chapter  4 

Finite  element  analysis 


4.1  Introduction 

The  use  of  finite  element  analysis  (FEA)  to  solve  problems  in  engineering  has  become  quite 
common  with  the  advent  of  many  commercially  available  programs.  The  magnetostrictive 
phenomenon  is  not  modeled  with  any  of  these  programs,  but  there  is  a  strong  analogy  be¬ 
tween  piezoelectricity  and  piezomagnetism.  When  this  analogy  is  fully  exploited,  piezomag¬ 
net  ic  behavior  can  be  modeled  with  the  piezoelectric  analysis  available  in  the  commercially 
available  program  ABAQUS.  The  geometry  of  the  model  presented  here  is  axisymmetric  and 
consists  of  a  Terfenol-D  rod  embedded  in  a  polymer  cylinder. 

4.1.1  Analogy  between  piezomagnetism  and  piezoelectricity 

The  governing  equations  for  piezomagnetism,  as  discussed  in  Section  1.2.2,  are  analogous  to 
those  of  the  piezoelectric  effect.  The  important  quantities  in  piezomagnetism  are  magnetic 
flux  (B)  and  magnetic  field  (H).  These  are  analogous  to  the  relevant  quantities  of  piezo¬ 
electricity,  electric  displacement  (D)  and  electric  field  (E),  respectively.  Table  4.1  shows 
a  comparison  of  the  governing  equations  for  piezomagnetism  and  piezoelectricity  [22],  To 
complete  the  analogy,  Table  4.2  shows  a  comparison  of  all  the  pertinent  quantities  and 
corresponding  units  [23,  22], 
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Table  4.1:  Comparison  of  piezomagnetic  and  piezoelectric  governing  equations. 
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Table  4.3:  Comparison  of  interface  conditions  for  piezomagnetism  and  piezoelectricity. 


Piezomagnetism 

Piezoelectricity 

Continuity  of  displacement 

u\  = 

u\  =  ulL 

Continuity  of  normal  stresses 

Continuity  of  scalar  potential 

0m  =  0m 

Continuity  of  normal  magnetic  flux 
or  normal  electric  displacement 

riiB\  =  riiBf 

n,D\  =  ntDf 

4.1.2  Boundary  conditions 

Along  with  the  piezomagnetic  governing  equations,  boundary  conditions  are  necessary  to 
solve  a  problem.  These  conditions  stem  from  the  equations  [10] 


V  •  B  =  0  (4.1a) 

V  x  H  =  0  (4.1b) 


which  hold  for  magnetostatic  problems  when  there  are  no  free  currents.  The  electrostatic 
equations 


V  •  D  =  0  (4.2a) 

V  x  E  =  0  (4.2b) 

hold  when  there  are  no  free  charges  [10].  These  analogous  sets  of  equations  lead  to  impor¬ 
tant  continuity  conditions  at  the  interface  of  two  materials.  A  comparison  of  the  interface 
conditions  for  piezomagnetism  and  piezoelectricity  follows  in  Table  4.3.  The  superscripts 
I  and  II  specify  the  two  materials  on  either  side  of  the  interface.  The  symbol  nt  denotes 
the  components  of  the  unit  normal  to  the  interface  and  rq  denotes  the  components  of  the 
displacement  vector. 

One  side  of  the  equation  for  continuity  of  normal  electric  displacement  is  often  set  to 
zero.  This  condition  is  called  the  air-dielectric  interface  condition  and  should  be  used  when 
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the  permittivity  of  one  material  is  much  larger  than  the  other  [23].  The  condition  could 
similarly  be  applied  to  the  magnetic  flux  for  appropriate  permeability  values.  Although  this 
approximation  is  a  common  one,  it  will  not  be  used  in  this  analysis. 

The  final  boundary  condition  is  setting  the  scalar  potential  equal  to  a  reference  potential 
on  the  surfaces  where  electrodes  are  located.  Although  these  clectroded  surfaces  do  not 
physically  exist  in  the  piezomagnetic  case,  the  locations  are  analogous  to  those  where  the 
magnetic  potential  is  set,  since  the  scalar  potentials  are  similarly  defined: 

V  (pe  =  E 
-V0m  =  H 

4.1.3  Symmetry  from  polarization 

When  Terfenol-D  is  not  polarized,  it  is  isotropic,  but  when  it  is  magnetized  (polarized), 
it  develops  transverse  isotropy  about  the  induced  polarization  vector  [32],  This  type  of 
symmetry  is  analogous  to  that  of  piezoelectric  polycrystallinc  ceramics  and  is  classified  as 
crystallographic  class  6mm  in  the  hexagonal  crystal  system.  The  induced  polarization  vector 
is  typically  denoted  by  A"3.  Any  direction  perpendicular  to  X3  can  be  defined  as  Xl}  since 
there  is  isotropy  in  the  plane  perpendicular  to  X3. 

The  presence  of  this  transverse  isotropy  about  the  3-axis  greatly  reduces  the  number  of 
independent  material  constants  necessary  to  characterize  Terfenol-D  fully.  The  form  of  the 


(4.3a) 

(4.3b) 
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clastic  compliance  matrix  becomes 


Sl  1  S13  S13  0  0  0 

512  Sn  S13  0  0  0 

513  Sl3  S33  0  0  0 

0  0  0  s44  0  0 

0  0  0  0  s44  0 

0  0  0  0  0  2(sn  —  S12) 

and  only  five  material  constants  are  needed  to  populate  it.  It  is  important  to  note  that  the 
compressed  matrix  form  has  been  used  and  the  superscripts  on  have  been  omitted.  The 
transverse  isotropy  also  modifies  the  piezomagnetic  matrix  to  the  form 

0  0  0  6?14  G?15  0 

0  0  0  di5  — rfu  0 

d'M  (h\  ^33  0  0  0 

with  four  material  constants.  Finally,  two  material  constants  fully  characterize  the  trans¬ 
versely  isotropic  permeability  matrix 

[An  0  0 

0  [in  0 
0  0  /r33 

Again,  superscripts  on  /i.^-  have  been  omitted  since  the  form  of  the  matrix  is  the  same  for 
both  cases.  In  summary,  the  matrix  form  of  the  piezomagnetic  equation  with  transverse 
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isotropy  about  X3  is 


Si 

oH 

*11 

sH 

*12 

*13 

0 

0 

0 

0 

0 

1 

CO 

^3 

T\ 

S2 

SH 
*  12 

sH 

*11 

SH 

*13 

0 

0 

0 

0 

0 

<^31 

t2 

S3 

SH 

*13 

SH 

*13 

*33 

0 

0 

0 

0 

0 

d33 

t3 

S4 

0 

0 

0 

sH 

*44 

0 

0 

di4 

di5 

0 

T4 

s5 

= 

0 

0 

0 

0 

*44 

0 

di5 

—d\4 

0 

T5 

Sq 

0 

0 

0 

0 

0 

2(*fi  -  *f2) 

0 

0 

0 

t6 

Bi 

0 

0 

0 

d  14 

dw 

0 

T 

/% 

0 

0 

Hi 

B-2 

0 

0 

0 

dis 

— di4 

0 

0 

T 

R11 

0 

h2 

b3 

^31 

d3i 

^33 

0 

0 

0 

0 

0 

T 

h33 

h3 

(4.4) 


4.2  Finite  element  formulation 


4.2.1  Preprocessing  with  PATRAN 

The  preprocessing  for  this  analysis  was  performed  using  PATRAN  6.2.  In  this  software,  the 
geometry  is  defined  first,  and  then  several  options  for  mesh  development  are  available.  The 
user  can  define  mesh  seeds  with  a  uniform  distribution,  a  one-way  bias,  or  a  two-way  bias 
along  curves.  The  program  then  generates  elements  and  nodes  from  these  mesh  seeds.  Two 
options,  isomesh  and  paver,  are  available  for  the  generation  of  the  finite  element  mesh  from 
the  mesh  seeds.  Isomesh  is  good  for  small  areas  where  the  user  wants  a  very  even  distribution 
of  elements.  Paver  works  wells  for  a  large  area  that  does  not  need  a  uniform  distribution  of 
the  elements.  In  general,  paver  generates  a  much  smaller  total  number  of  elements. 

After  the  mesh  is  developed,  the  user  must  input  the  boundary  conditions,  loading  con¬ 
ditions,  and  material  properties.  PATRAN  6.2  does  not  have  a  degree  of  freedom  for  electric 
potential  available,  so  the  analysis  can  not  be  completed  with  this  program.  The  final  step 
is  to  produce  an  input  file  for  ABAQUS. 
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Table  4.4:  Piezoelectric  elements  available  in  ABAQUS. 


Plane  strain 

elements 

Plane  stress 

elements 

Axisymmetric 

elements 

Three-dimensional 

elements 

CPE3E 

CPS3E 

CAX3E 

C3D4E 

CPE4E 

CPS4E 

CAX4E 

C3D6E 

CPE6E 

CPS6E 

CAX6E 

C3D8E 

CPE8E 

CPS8E 

CAX8E 

C3D10E 

CPE8RE 

CPS8RE 

CAX8RE 

C3D15E 

C3D20E 

C3D25E 

4.2.2  Analysis  with  ABAQUS 

ABAQUS,  Version  5.7,  has  piezoelectric  elements  and  can  perform  piezoelectric  analyses. 
Procedures  are  available  for  both  static  and  dynamic  piezoelectric  analyses.  The  piezoelec¬ 
tric  analysis  procedure  uses  the  equations  presented  in  Table  4.1.  Since  these  equations  and 
the  boundary  conditions  in  Table  4.3  have  a  one-to-one  correspondence  with  those  for  piezo- 
magnetism,  one  can  effectively  model  piezomagnetism  by  substitution  of  the  appropriate 
quantities  into  the  piezoelectric  model  via  Table  4.2. 

There  are  several  types  of  piezoelectric  elements  available  in  ABAQUS.  Plane  strain, 
plane  stress,  and  axisymmetric  elements  are  available  with  3,  4,  6,  and  8  nodes.  Three- 
dimensional  elements  are  available  with  4,  6,  8,  10,  15,  and  20  nodes.  The  piezoelectric 
elements  have  degrees  of  freedom  for  displacement  and  scalar  electric  potential.  The  elements 
are  named  according  to  a  certain  pattern.  The  first  letter  corresponds  to  the  major  type  of 
element.  For  piezoelectric  analysis,  the  first  letter  is  always  C ,  which  stands  for  continuum. 
The  next  two  characters  correspond  to  plane  strain  (PE),  plane  stress  (PS),  axisymmetric 
(AX’),  or  three-dimensional  (3D).  The  next  number  is  the  number  of  nodes  per  element. 
After  this  number,  such  options  as  piezoelectric  (E)  and  reduced  integration  (R)  are  included. 
Table  4.4  shows  a  complete  listing  of  the  piezoelectric  elements  available  in  ABAQUS. 

Boundary  conditions  can  be  applied  to  displacement  and  electric  potential  at  any  nodes. 
Displacement  degrees  of  freedom,  ut,  are  numbered  1,  2,  and  3.  Electric  potential,  0e,  is 
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numbered  degree  of  freedom  9.  Loading  is  applied  to  an  element  or  an  element  face.  Body 
force,  edge  pressure,  surface  pressure,  and  electric  charge  are  several  of  the  loading  quantities 
available  for  a  piezoelectric  analysis. 

The  outputs  from  the  piezoelectric  analysis  include  both  nodal  quantities  and  integration 
point  quantities.  The  nodal  quantities  are  electric  potential  and  displacement.  Integration 
point  values  are  not  calculated  at  the  nodes  of  an  element,  but  rather  at  integration  points 
inside  the  boundaries  of  an  element.  Integration  point  variables  include  stress,  strain,  electric 
displacement,  and  electric  held. 

4.2.3  Modification  of  the  input  file 

After  the  PATRAN  input  hie  is  created,  it  must  be  modified  before  it  is  ready  for  piezoelectric 
analysis  in  ABAQUS.  The  input  and  analysis  are  treated  as  piezoelectric  by  ABAQUS,  but 
if  the  user  enters  the  appropriate  quantities  for  piezomagnetism,  the  results  will  correspond 
to  piezo  magnetism.  The  analogy  of  relevant  quantities  is  completely  described  in  Table  4.2. 

An  abbreviated  example  input  hie  is  included  in  Appendix  C.  The  hrst  item  to  modify 
in  the  PATRAN  hie  is  the  type  of  element.  In  the  example  hie,  PATRAN  created  CAX8 
elements.  For  piezoelectric  analysis,  the  element  type  must  be  changed  to  CAX8E,  which 
denotes  axisymmetric,  8-node,  piezoelectric  elements. 

Second,  the  material  properties  must  be  adjusted  to  include  piezomagnetic  properties. 
PATRAN  produces  the  input  only  for  materials  with  elastic  properties,  which  are  presented 
in  the  order  of  Young’s  modulus  followed  by  Poisson’s  ratio.  The  piezomagnetic  coefficients 
can  be  entered  as  the  d  coefficients  or  the  e  coefficients  described  in  Table  4.2.  If  the  user 
enters  d  coefficients,  “type  =  e”  must  be  included.  For  e  coefficients,  “type  =  s”  is  the 
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necessary  flag.  The  format  of  the  entry  of  the  coefficients  is 


dm, 

dl22, 

dl33, 

dll2, 

^113, 

to 

CO 

^211,  C?222 

^233, 

^212, 

^213, 

d 223, 

^311, 

dz22 , 

^333,  ^312 

^313,  ^323 

and  this  pattern  must  be  followed  exactly.  It  is  important  to  note  that  this  order  is  slightly 
different  from  the  order  of  the  coefficients  in  the  matrices  described  in  the  IEEE  Stan¬ 
dards  [22,  23]  and  in  Equation  (1.2). 

The  permittivity  coefficients  of  piezoelectricity  are  referred  to  as  dielectric  properties 
in  ABAQUS.  There  is  no  flag  available  for  constant  stress  or  constant  strain,  so  the  user 
may  input  whichever  properties  are  appropriate  to  the  analysis.  For  piezomagnetism,  the 
analogous  permeability  quantities  are  input.  They  can  be  entered  in  the  order 

Pll,  P-12,  P 22,  pl3,  p23,  p33 

with  “type  =  aniso” .  If  symmetry  exists,  “type  =  ortho”  and  “type  =  iso”  are  possible  flags 
that  reduce  the  number  of  quantities  that  need  to  be  included.  The  permeability  values 
should  not  be  entered  as  relative  permeability  with  respect  to  po  but  must  be  entered  as 
true  permeability. 

PATRAN  can  not  produce  boundary  conditions  for  the  scalar  potential,  so  the  input  file 
must  be  modified  to  include  these.  Electric  flux  loads  must  also  be  input  by  hand  since  they 
do  not  exist  in  the  initial  PATRAN  file.  Finally,  the  line RESTART.  WRITE”  must  be 
added  to  the  end  of  the  input  file  just  before  “END  STEP”.  This  command  creates  a  file 
from  which  ABAQUS  Post  can  run. 
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4.3  Composite  cylinder  model 


4.3.1  Geometry,  loading,  and  boundary  conditions 

The  geometry  of  concentric  cylinders  was  selected  to  model  a  Terfenol-D  rod  embedded  in 
a  polymer  matrix.  As  shown  in  Figure  4.1,  the  piezomagnetic  rod  has  a  radius  that  is  one- 
tenth  of  the  outer  radius  of  the  polymer  cylinder,  giving  a  volume  fraction  of  1%  for  the 
rod.  The  rod’s  poling  axis  is  parallel  to  the  axis  of  the  cylinder  (2-direction).  The  loading 
is  a  uniform  compressive  stress  of  —10  MPa  in  the  2-direction  applied  to  the  ends  of  the 
cylinder.  Although  this  geometry  is  three-dimensional,  it  can  be  modeled  as  axisymmetric 
by  choosing  the  boundary  conditions  appropriately. 

The  left  edge  of  the  model  is  a  symmetry  boundary,  so  the  normal  displacement  is 
constrained  and  U\  (displacement  in  the  1-direction)  is  set  equal  to  zero  along  this  edge.  In 
the  piezoelectric  case,  the  electric  potential  would  be  set  to  zero  along  the  top  and  bottom 
of  the  cylinder.  The  analogous  piezomagnetic  boundary  condition  is  to  set  the  magnetic 
potential,  0m,  equal  to  zero  on  top  and  bottom  edges.  The  boundary  conditions  are  also 
denoted  in  Figure  4.1. 

4.3.2  Material  properties 

Several  assumptions  are  made  about  the  properties  of  Terfenol-D  in  order  to  complete  the 
model.  First,  since  only  a  few  clastic  properties  of  Terfenol-D  are  published,  the  assumption 
of  isotropic  elasticity  is  assumed.  A  value  of  50  GPa  is  used  for  Young’s  modulus  [32]  and 
0.444  is  used  for  the  Poisson’s  ratio  [19].  For  the  piezomagnetic  coefficients,  the  “type  = 
e”  flag  is  used.  The  value  of  d 33  is  taken  to  be  1.48  x  10-10  m/A,  using  the  results  from 
Section  2.4.3. 1  and  assuming  transverse  isotropy  with  respect  to  the  3-axis.  Three  more 
piezomagnetic  coefficients  are  necessary  to  populate  the  piezomagnetic  matrix.  The  values 
for  d3i  and  d35  are  assumed  to  be  0.5d33  and  — 0.5d33,  respectively.  When  a  magnetostrictive 
material  is  poled  in  the  3-direction,  du  has  not  been  experimentally  observed  [32],  so  it 


Table  4.5:  Summary  of  material  properties  used  in  finite  element  analysis  of  composite 
cylinder  model. 


E  (GPa) 

V 

d33  (m/A) 

d3 1  (m/A) 

di5  (m/A) 

pa  (H/m) 

Polyester 

3.8 

0.33 

1.48  x  10“12 

-0.72  x  10”12 

0.72  x  10“12 

6.28  x  10“6 

Terfenol-D 

50 

0.444 

0 

0 

0 

1.26  x  10”6 

is  taken  to  be  zero.  ABAQUS  defines  the  radial  direction  as  the  1 -direction  and  the  axial 
direction  as  the  2-direction  for  axisymmetric  elements.  For  this  model,  the  axial  direction 
is  the  poling  axis,  so  the  piezomagnetic  coefficients  must  be  transformed  to  the  coordinate 
system  of  the  model  shown  in  Figure  4.1. 

4.4  Convergence  study 

Two  different  meshes  were  constructed  for  the  finite  element  analysis  of  the  composite  cylin¬ 
der.  The  elements  used  were  axisymmetric  8-node  quads  (CAX8E)  for  both  meshes.  Each 
node  has  two  displacement  degrees  of  freedom,  u\  and  u2 ,  and  a  degree  of  freedom  for  electric 
(in  this  case,  magnetic)  potential,  (j)e.  Mesh  1,  shown  in  Figure  4.2,  has  14,869  nodes  and 
4,  860  elements.  This  mesh  was  created  using  the  isomesh  option  in  PATRAN.  A  second 
mesh  was  constructed  with  a  higher  concentration  of  small  elements  at  the  interface.  Mesh 
2,  shown  in  Figure  4.3,  has  17,  750  nodes  and  5,  813  elements.  The  mesh  in  the  Terfenol-D 
rod  was  created  with  isomesh  and  the  mesh  in  the  polyester  matrix  was  created  with  paver. 

To  determine  how  well  each  of  these  meshes  modeled  the  composite  cylinder  problem, 
both  meshes  were  run  and  the  results  were  compared.  The  shear  stress,  aq2,  and  the  magnetic 
flux,  Bi,  were  measured  at  the  interface  of  the  Terfenol-D  rod  and  the  polyester  matrix. 
These  two  quantities  should  be  continuous  across  the  interface  as  shown  in  Table  4.3.  Using 
ABAQUS  Post,  the  stress  and  flux  could  be  determined  for  each  side  of  the  interface.  At 
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the  interface,  the  equations 


Dj  -  Bl  =  0  (4.5a) 

O12  -  a12  =  0  (4.5b) 

should  be  true,  where  the  superscript  T  denotes  the  Terfenol-D  side  of  the  interface  and  P 
denotes  the  polyester  side  of  the  interface.  For  both  meshes,  the  continuity  was  enforced 
except  at  the  top  and  bottom  of  the  model.  Figures  4.4  and  4.5  show  the  difference  in  B\ 
across  interface  for  the  lower  and  upper  10%  of  the  interface,  respectively.  Both  meshes 
enforce  the  continuity  well  up  to  approximately  5%,  where  mesh  1  begins  to  fail.  Mesh  2  is 
able  to  capture  the  continuity  much  better.  Only  1%  of  the  interface  at  the  top  and  bottom 
of  the  model  has  problems  with  continuity.  Similar  results  are  seen  for  the  continuity  of 
shear  stress  in  Figures  4.6  and  4.7.  The  increased  number  and  decreased  size  of  the  elements 
at  the  interface  in  Mesh  2  significantly  improved  the  model  results. 

4.5  Results 

All  results  shown  were  obtained  using  Mesh  2  unless  otherwise  specified.  The  variables  of 
interest  are  the  two  components  of  magnetic  flux.  ABAQUS  Post  was  used  to  obtain  contour 
plots  of  the  magnetic  flux.  The  model  was  separated  into  two  parts — the  Terfenol-D  rod  and 
the  polyester  matrix,  for  the  contour  plots.  Figures  4.8  and  4.9  show  the  magnetic  fluxes,  B\ 
and  B2,  in  the  Terfenol-D  rod.  The  polyester  matrix  contour  plots  are  given  in  Figures  4.10 
and  4.11.  The  values  for  the  contours  were  adjusted  for  each  plot  so  that  more  detail  could 
be  seen. 

The  magnetic  flux,  B\ ,  is  continuous  across  the  interface,  so  the  value  remains  rela¬ 
tively  large  in  the  matrix  near  the  interface,  but  drops  off  quickly  as  the  distance  from  the 
Terfenol-D  rod  increases.  Only  2  mm  from  the  rod,  the  flux  is  nearly  zero.  The  magnetic 


flux,  B2,  is  not  continuous  across  the  interface.  In  the  Terfenol-D  rod,  B2  is  larger  than 
any  other  magnetic  flux  in  the  model.  At  the  midpoint  of  the  rod  near  the  interface,  B2  is 
around  —72  G.  Across  the  interface  in  the  matrix,  B2  is  an  order  of  magnitude  smaller  and 
of  opposite  sign.  The  quick  decay  of  the  magnetic  flux  in  the  matrix  is  demonstrated  clearly 
in  Figures  4.12  and  4.13.  These  two  figures  show  the  magnetic  flux  along  the  line  y  =  5  mm. 

To  compare  the  finite  element  model  results  with  the  compression  test  results  presented 
in  Section  2.4.1,  the  magnetic  flux  along  the  interface  in  the  Terfenol-D  rod  is  plotted  in  Fig¬ 
ures  4.14  and  4.15.  The  trends  are  similar  to  the  gage  length  results  for  the  pure  Terfenol-D 
specimen.  The  “fanning  phenomenon”  is  seen  in  the  magnetic  flux,  Bi}  that  is  perpendicular 
to  the  direction  of  the  magnetization  or  poling  axis.  Also,  the  response,  B 1;  in  the  direction 
of  poling  is  nearly  constant  in  the  center  of  the  specimen,  but  decreases  in  magnitude  slightly 
before  rising  dramatically  at  the  ends  of  the  specimen.  Both  these  trends  are  documented 
in  Figure  2.13. 

There  is  a  large  difference  in  the  magnitude  of  the  magnetic  flux  measured  in  the  com¬ 
pression  testing  and  the  magnetic  flux  calculated  in  the  finite  element  model.  The  large 
discrepancy  is  seen  mainly  in  the  magnetic  flux  along  the  loading  axis,  which  is  also  the 
poling  axis.  The  finite  element  model  uses  linear  piezomagnetism,  which  is  not  a  poor 
assumption,  but  there  is  no  concept  of  saturation  incorporated  into  the  model.  If  the  piezo- 
magnetic  analysis  were  modified  to  include  saturation  behavior,  the  exaggerated  predictions 
of  magnetic  flux  would  disappear. 
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Figure  4.1:  Geometry  of  the  composite  cylinder  model. 


90 


2 

A 


91 


Figure  4.3:  Mesh  2  with  17,  750  nodes  and  5,  813  elements. 
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Figure  4.4:  Difference  in  magnetic  flux,  B i,  across  interface  at  lower  10%  of  the  composite 
cylinder  model. 


Position  along  interface,  y  (mm) 


Figure  4.5:  Difference  in  magnetic  flux,  B\ ,  across  interface  at  upper  10%  of  the  composite 
cylinder  model. 
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Figure  4.6:  Difference  in  shear  stress,  <712,  across  interface  at  lower  10%  of  the  composite 
cylinder  model. 


Figure  4.7:  Difference  in  shear  stress,  012,  across  interface  at  upper  10%  of  the  composite 
cylinder  model. 
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Figure  4.8:  Contour  plot  of  the  magnetic  flux,  B i,  in  the  Terfenol-D  rod  of  the  composite 
cylinder  model. 
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Figure  4.9:  Contour  plot  of  the  magnetic  flux,  B 2,  in  the  Terfenol-D  rod  of  the  composite 
cylinder  model. 


96 


2 

A 


Bj(G) 

-3.21 

-1.60 

-1.31 

-1.02 

-0.73 

-0.44 

-0.15 

+0.15 

+0.44 

+0.73 

+1.02 

+1.31 

+1.60 

+3.25 


Figure  4.10:  Contour  plot  of  the  magnetic  flux,  B i,  in  the  polyester  matrix  of  the  composite 
cylinder  model. 
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Figure  4.11:  Contour  plot  of  the  magnetic  flux,  B 2,  in  the  polyester  matrix  of  the  composite 
cylinder  model. 


Figure  4.12:  The  magnetic  flux,  B i,  along  the  line  y  =  5  mm  of  the  composite  cylinder 
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Figure  4.13:  The  magnetic  flux,  B2,  along  the  line  y 
model. 
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Figure  4.14:  Magnetic  flux,  B f,  along  the  interface  in  the  Terfenol-D  rod. 
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Magnetic  flux,  B%,  along  the  interface  in  the  Terfenol-D  rod. 
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Chapter  5 

Summary  and  conclusions 


The  increasing  use  of  composite  materials  in  structural  applications  has  led  to  a  need  to 
develop  new  structural  health  monitoring  techniques.  Tagging  of  a  structure  with  a  smart 
material  phase  is  a  viable  option.  The  work  reported  in  this  thesis  was  undertaken  to 
understand  the  magnetomechanical  behavior  of  Terfenol-D  and  Terfenol-D  particulate  com¬ 
posite  to  further  the  possibility  of  composites  tagged  with  the  magnetostrictive  material, 
Terfenol-D. 

Uniaxial  compression  testing  was  conducted  on  a  pure  Terfenol-D  specimen  and  several 
Terfenol-D  particulate  composites.  The  composites  consisted  of  a  polyester  matrix  and  2.2% 
volume  fraction  of  Terfenol-D  particulates.  Different  combinations  of  curing  magnetic  flux 
and  annealing  magnetic  flux  were  studied.  Chaining  of  the  particles  was  seen  whether  the 
specimen  was  cured  with  a  flux  through  its  length  or  its  width.  The  three  components  of 
magnetic  flux  were  measured  with  a  Hall  probe  over  the  gage  length  of  each  specimen  under 
uniaxial  compressive  loading.  Results  showed  that  maximum  responses  occurred  when  the 
maximum  number  of  domains  were  oriented  in  the  direction  of  the  applied  compressive  stress. 
It  was  shown  that  if  the  domains  are  initially  aligned  in  the  direction  of  the  compressive 
stress,  they  rotate  away  from  this  axis  as  the  stress  increases.  Also,  a  “fanning  phenomenon” 
was  observed  that  could  be  explained  by  observing  the  magnetic  field  lines  of  a  uniformly 
magnetized  cylinder. 

Piezomagnetic  coefficients  were  measured  for  the  Terfenol-D  specimen  and  the  compos- 
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ites.  The  largest  piezomagnetic  coefficient  measured  was  d33  (1.48  x  1CT10  m/A)  for  the  pure 
specimen.  This  value  is  much  lower  than  the  values  reported  in  the  literature  for  the  direct 
effect.  Saturation  behavior  was  observed  for  the  pure  specimen  at  stresses  as  low  as  —2  MPa. 
The  values  of  the  piezomagnetic  coefficients  were  at  least  two  orders  of  magnitude  smaller 
for  the  Terfenol-D  composites.  The  largest  d33  coefficient  measured  for  a  Terfenol-D  compos¬ 
ite  was  7.86  x  10~12  m/A.  Saturation  was  not  always  reached  in  the  Terfenol-D  composites 
tested. 

Cyclic  testing  was  conducted  on  the  Terfenol-D  composites.  First,  a  specimen  was  mag¬ 
netically  annealed.  Then,  it  was  loaded  and  unloaded  in  compression.  The  loading  and 
unloading  cycle  was  repeated  until  a  stable  response  was  seen  without  any  further  anneal¬ 
ing.  It  was  shown  that  as  few  as  3  cycles  were  necessary  to  reach  a  stable  response.  The  cyclic- 
testing  showed  that  increasing  the  magnitude  of  the  curing  magnetic  flux  increased  both  the 
initial  and  the  stable  responses.  On  the  other  hand,  increasing  the  annealing  magnetic  flux 
increased  only  the  initial  magnetic  response. 

A  possible  explanation  for  the  absence  of  the  saturation  behavior  in  the  composites  is 
that  the  residual  stresses  are  quite  high.  To  investigate  the  residual  stress  state,  a  modified 
Eshelby  model  was  developed.  Thermal  mismatch,  chemical  shrinkage,  and  magnetostrictive 
strains  were  all  included  in  the  Eshelby  model.  Two  different  geometries  were  studied: 
a  single  particle  modeled  as  a  spherical  inclusion  and  a  chain  of  particles  modeled  as  an 
ellipsoidal  inclusion. 

A  finite  element  analysis  was  conducted  on  a  composite  cylinder  in  compression  using  the 
piezoelectric  analysis  available  in  ABAQUS.  Results  from  the  finite  element  analysis  showed 
the  same  “fanning  phenomenon”  seen  in  the  compression  tests.  Also,  the  magnetic  flux 
change  was  largest  in  the  direction  of  the  poling  axis,  just  as  in  the  compression  testing.  The 
limitation  to  using  the  finite  element  analysis  is  that  saturation  behavior  is  not  currently 
incorporated  into  the  governing  equations.  The  compression  testing  of  the  pure  Terfenol-D 
rod  showed  that  saturation  occurred  at  very  low  stresses,  so  the  finite  element  analysis  greatly 
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overpredicted  the  magnetic  response  observed  in  experiments. 

This  research  could  be  improved  with  a  better  magnetic  flux  sensing  device.  Such  a  device 
would  simultaneously  measure  all  three  components  of  the  magnetic  flux  at  the  same  point. 
Unfortunately,  the  measurements  of  flux  in  this  research  were  not  all  taken  at  the  same 
distance  from  the  surface  of  the  specimen.  It  was  shown  that  the  magnetic  flux  measured  is 
strongly  affected  by  the  distance  from  the  surface  of  the  specimen.  Calculations  of  maximum 
magnetic  flux  change  and  piezomagnetic  coefficients  should  be  compared  carefully  as  some  of 
the  measurements  are  smaller  since  they  were  taken  farther  from  the  surface  of  the  specimen. 

The  finite  element  analysis  could  be  improved  by  incorporating  saturation  behavior  into 
the  governing  equations  used  in  the  piezoelectric  analysis.  Also,  it  would  be  interesting  to 
include  the  residual  stresses  predicted  by  the  modified  Eshelby  model.  These  two  modifi¬ 
cations  would  probably  improve  the  correlation  between  the  experimental  results  and  the 
finite  element  analysis  results. 
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Appendix  A 

Gage  length  characterization  results 
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Figure  A.l:  Magnetic  flux  ,B i,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 
flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 


Compressive  stress,  on  (MPa) 


Figure  A. 2:  Zeroed  magnetic  flux,  B i,  over  the  gage  length  of  pure  Terfenol-D  specimen, 
annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 3:  Magnetic  flux,  B 2,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 
flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 


Compressive  stress,  on  (MPa) 


Figure  A. 4:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  pure  Terfenol-D  specimen, 
annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 5:  Magnetic  flux,  B 2,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 
flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 


Compressive  stress,  on  (MPa) 


Figure  A. 6:  Zeroed  magnetic  flux,  B2 ,  over  the  gage  length  of  pure  Terfenol-D  specimen, 
annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 
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Figure  A. 7:  Magnetic  flux,  B%,  over  the  gage  length  of  pure  Terfenol-D  specimen,  annealing 
flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 


Compressive  stress,  on  (MPa) 


Figure  A. 8:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  pure  Terfenol-D  specimen, 
annealing  flux  of  —800  G  in  1-direction  applied  before  testing  at  each  position. 
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1.3 


Compressive  stress,  on  (MPa) 


Figure  A. 9:  Magnetic  flux,  B i,  over  the  gage  length  of  (1,1,1)  Tcrfenol-D  composite  specimen, 
curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each  position. 


Compressive  stress,  on  (MPa) 


Figure  A. 10:  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 11:  Magnetic  flux,  B2 ,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 
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Compressive  stress,  an  (MPa) 


Figure  A. 12:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 13:  Magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,1,1)  Terfenol-D 
composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position. 


Compressive  stress,  an  (MPa) 


Figure  A. 14:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,1,1) 
Terfenol-D  composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied 
before  testing  at  each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 15:  Magnetic  flux,  B3,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 


Figure  A. 16:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (1,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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0.5 


Compressive  stress,  on  (MPa) 


Figure  A. 17:  Magnetic  flux,  B i,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 


Figure  A. 18:  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 19:  Magnetic  flux,  B2 ,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Figure  A. 20:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 21:  Magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,3,1)  Tcrfenol-D 
composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before 
testing  at  each  position. 


Compressive  stress,  an  (MPa) 


Figure  A. 22:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  on  the  side  face  of  (1,3,1) 
Terfenol-D  composite  specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied 
before  testing  at  each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 23:  Magnetic  flux,  B3,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 


Figure  A. 24:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (1,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 

Figure  A. 25:  Magnetic  flux,  B i,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 
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Compressive  stress,  on  (MPa) 

Figure  A. 26:  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 27:  Magnetic  flux,  B2)  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  Gn  (MPa) 


Figure  A. 28:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 

Figure  A. 29:  Magnetic  flux,  B3,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 30:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (3,1,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 31:  Magnetic  flux,  B i,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Figure  A. 32:  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 

Figure  A. 33:  Magnetic  flux,  B2)  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 

Figure  A. 34:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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0.05 


Compressive  stress,  on  (MPa) 


Figure  A. 35:  Magnetic  flux,  B3,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Figure  A. 36:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (3,2,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 37:  Magnetic  flux,  B i,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 


Figure  A. 38:  Zeroed  magnetic  flux,  Bi,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 39:  Magnetic  flux,  B2)  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Figure  A. 40:  Zeroed  magnetic  flux,  B2,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Compressive  stress,  on  (MPa) 


Figure  A. 41:  Magnetic  flux,  B3,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite  spec¬ 
imen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at  each 
position. 


Compressive  stress,  on  (MPa) 


Figure  A. 42:  Zeroed  magnetic  flux,  B3,  over  the  gage  length  of  (3,3,1)  Terfenol-D  composite 
specimen,  curing  flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  testing  at 
each  position. 
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Appendix  B 
Cyclic  testing  results 
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Figure  B.l:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 


Figure  B.2:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 
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Figure  B.3:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 


Figure  B.4:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 
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-2.4 


Compressive  stress,  on  (MPa) 


Figure  B.5:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 


Figure  B.6:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 


131 


Figure  B.7:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 


Figure  B.8:  Cyclic  magnetic  response,  B i,  of  (1,1,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 
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Figure  B.9:  Cyclic  magnetic  response,  B3,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 


Figure  B.10:  Cyclic  magnetic  response,  B%,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —800  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 
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Figure  B.ll:  Cyclic  magnetic  response,  B%,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 


Figure  B.12:  Cyclic  magnetic  response,  B%,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —800  G  applied  before  first  cycle  only. 
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Compressive  stress,  on  (MPa) 


Figure  B.13:  Cyclic  magnetic  response,  B%,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 


Figure  B.14:  Cyclic  magnetic  response,  B%,  of  (3,3,1)  Terfenol-D  composite  specimen,  curing 
flux  of  —2600  G,  and  annealing  flux  of  —2600  G  applied  before  first  cycle  only. 
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Appendix  C 

Sample  input  file  for  FEA 


❖HEADING,  SPARSE 
** 

❖NODE 

1,  0.000102848,  0. 

2,  0.0001,  0. 

(remaining  node  list  omitted) 

** 

** 

❖ELEMENT,  TYPE=CAX8E,  ELSET=R0D 

3601,  14789,  14626,  14621,  14784, 

14785,  14791 

3602,  14626,  14462,  14457,  14621, 

14619,  14627 

(remaining  element  list  omitted) 

** 

** 

** 

❖ELEMENT,  TYPE=CAX8E,  ELSET=MATRIX 

1,  2,  3,  167,  166, 

165,  4 

2,  166,  167,  331,  330, 

329,  274 

(remaining  element  list  omitted) 

** 

** 

** 

❖❖  chain 
** 

❖SOLID  SECTION, 

ELSET=CHAIN,  MATERI AL=TERFEN0L 

1., 

** 

❖❖  matrix 
** 

❖SOLID  SECTION,  ELSET=MATRIX ,  MATERI AL=P0LYMER 

1., 

** 

** 


14790, 

14625, 


1, 

165, 


14627, 

14463, 


5, 

275, 
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❖MATERIAL,  NAME=POLYMER 
** 

❖ELASTIC,  TYPE=ISO 

3 . 8E+9 ,  0.33 

** 

❖PIEZOELECTRIC ,  TYPE=E 

0. ,0. ,0. ,0. ,0. ,0. ,0. ,0. , 

0. ,0. ,0. ,0. ,0. ,0. ,0. ,0. , 

0.  ,0. 

** 

❖DIELECTRIC,  TYPE=ISO 
12.56E-7 
** 

** 

** 

❖MATERIAL,  NAME=TERFENOL 
** 

❖ELASTIC,  TYPE=ISO 

5.E+10,  0.444 

** 

** 

❖PIEZOELECTRIC,  TYPE=E 

0. ,0. ,0. ,0.74E-10,0. ,0. , -0 . 74E-10 , 1 . 48E-10 , 

-0 . 74E-10 ,0.,0.,0.,0.,0.,0.,0., 

0. .0.74E-10 
** 

❖DIELECTRIC,  TYPE=ISO 

6.28E-6 

** 

** 

❖❖  noxdisp 
** 

❖BOUNDARY,  OP=NEW 

85,  1,,  0. 

220,  1,,  0. 

(remaining  nodes  and  boundary  conditions  omitted) 

** 

❖❖  potential 
** 

❖BOUNDARY,  0P=NEW 
POTENTIA,  9,,  0. 

** 

❖❖  Step  1,  Default  Static  Step 
❖❖  LoadCase,  Default 
** 

❖STEP,  AMPLITUDE=RAMP ,  PERTURBATION 
Linear  Static  Analysis 

❖❖  This  load  case  is  the  default  load  case  that  always  appears 
** 

❖STATIC 

** 

❖NSET,  NSET=POTENTIA 

1,  2,  3,  6,  7,  9,  10,  12, 
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(remaining  node  set  omitted) 

** 

** 

*ELSET,  ELSET=COMPRESS 

1,  91,  181,  271,  361,  451,  541, 

(remaining  element  set  omitted) 

** 

** 

**  compress 
** 

*DL0AD,  0P=NEW  COMPRESS, 

PI,  l.E+7 

** 

* RESTART,  WRITE 
** 

*END  STEP 


631, 
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